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Abstract 


This report describes measurements of the dielectric constant of seawater at a frequency of 1 .413 
GHz that is at the center of the L-Band radiometric protected frequency spectrum. Aquarius will 
be sensing the sea surface salinity from space in this band. The objective of the project is to 
refine the model function for the dielectric constant as a function of salinity and temperature so 
that remote sensing measurements can be made with the accuracy needed to meet the 
measurement goals (0.2 psu) of the Aquarius mission. 

The measurements were made, using a microwave cavity operated in the transmission 
configuration. The cavity’s temperature was accurately regulated to 0.02°C by immersing it in a 
temperature controlled bath of distilled water and ethanol glycol. Seawater had been purchased 
from Ocean Scientific International Limited (OSIL) at salinities of 30, 35 and 38 psu. 
Measurements of these seawater samples were then made over a range of temperatures, from 
10°C to 35°C in 5°C intervals. Repeated measurements were made at each temperature and 
salinity. Mean values and standard deviations were then computed. Total error budgets indicated 
that the real and imaginary parts of the dielectric constant had a relative accuracy of about 1%. 
Comparisons with the model function of Klein and Swift showed that the real part of the 
dielectric constant was about 1 .5% below their values while the imaginary part, on the average, 
was about 2 % below the Klein and Swift results. The differences in the imaginary parts between 
the two results increased with rising temperature and salinity. 
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1 Introduction 


This report describes L band measurements of the dielectric constant of seawater using a 
resonant cavity approach. The measurements have been made in support of the space-borne 
Aquarius mission with its goal of measuring sea surface salinity (SSS) to an accuracy of 0.2 psu. 
Because of the small variation in sea surface salinity over the open ocean, highly accurate L band 
radiometric measurements are required to sense these variations. Sensitivity calculations show 
that in order to achieve the required accuracy of 0.2 psu for sea surface salinity measurements, 
the brightness temperature measured by a satellite radiometer must be accurate to 0.3° K. This 
requirement, in turn, places an accuracy requirement of 0.2 % on the dielectric constant 
measurements used in the retrieval (i.e., the construction of the model function, e(T,S)). It is 
not clear that existing model functions such as Klein and Swift (1977) satisfy these accuracy 
requirements. The Klein and Swift model function, for example, is based, in part, on the 
dielectric measurement data of Ho and Hall (1973) and Ho et al (1974). These measurements 
were made over thirty years ago without the advantages of modern test equipment. In addition, 
they did not adequately cover the range of temperatures and salinity of interest to Aquarius 
mission. 

Salinity measurements made before 1978 were expressed in parts per thousand (%o) while after 
1978 the measurements were expressed in practical salinity units (psu). The model function of 
Klein and Swift was based on measurements made in parts per thousand. The new 
measurements, presented in this report, were made in practical salinity units. Comparison of the 
data presented in this report with the Klein and Swift model function will not be a problem, 
however. Lewis and Perkins (1981) point out that the two scales differ by less than .001 in the 
salinity range, 25<S<38 for salinities computed from chlorinity using the Knudsen relationship 
(Forch, et al 1902). Since the Ho et al (1974) measurements were made in parts per thousand 
using the Knudsen relationship, measurements made in parts per thousand and in practical 
salinity units can be used interchangeably with only negligible error. 

The dielectric constant measurements of seawater described in this report have been made at 
1.413 GHz. This is approximately the center frequency of the quiet zone (radiometric protected 
band) from 1 .400 GHz to 1 .427 GHz that Aquarius will utilize when making its measurements. 
The measurements have been made at temperatures of 10°C to 35°C in 5° increments with 
salinity values of 30, 35 and 38 psu. Although the range of the measurement set should be 
extended, these measurements represent values characteristic of open ocean conditions - the 
primary objective of the Aquarius mission. 

In this report a time dependence of exp(-zVuf) will be assumed. This assumed time dependence 
for Maxwell’s time harmonic equations leads to the expression of the complex relative dielectric 
constant as £-£ +ie " where £ and £*are the real and imaginary parts of the relative dielectric 
constant respectively. For lossless or lossy media, e" will be greater or equal to zero (£ r > 0). 

The report is structured as follows: in Section 2, sensitivity requirements for the dielectric 
constant measurements of seawater are derived based on the measurement accuracy goal of 0.2 
psu. In the same section, a review of past measurement data for this frequency band will be 


made. In addition, the use of out of band data to estimate in-band dielectric constants using 
model functions will be presented. In Section 3, a review of the available methods for the 
accurate measurement of the dielectric constant of lossy liquids is given, along with a detailed 
discussion of how past measurements have been made. In Section 4 the resonant cavity approach 
used by The George Washington University (GW) group is described along with the 
measurement data and the method used for calibration of the cavity. Finally in Section 5, 
conclusions and future work are discussed. 

2 Background 

In the following two subsections, there is presented a brief discussion of how the brightness 
temperature - measured over a half space of seawater - varies with temperature and salinity. In 
this context, the sensitivity of salinity to errors in the complex dielectric constant of seawater is 
then explored. This analysis is followed by a survey of seawater model function development 
over the last 40 years. 

2.1 Remote Sensing 

In remote sensing of sea surface salinity, the observed brightness temperature is directly related 
to the dielectric constant of the sea surface. The dielectric constant of the seawater is then related 
to the salinity, S, and the temperature, T, of the seawater by the model function, e(T,S ) . These 
relationships are illustrated in eq(l ) for the case of a radiometer looking directly down, normal to 
a flat sea surface (the simplest case). Here the brightness temperature, T B , is given by 

T„ = (l-|/?| 2 )r (la) 

where R is the Fresnel reflection of the surface, that is, 

R = (\-yfe)/(\+^). (lb) 

The relative dielectric constant of the seawater, e , is given in terms of the model function: 

£ = £(T,S) (lc) 

The dielectric constant contains the electrical properties of the medium and it is here that the 
dependence on salinity appears. 

As an illustration, eq(la) is shown in Figure 1 using the model function developed by Klein and 
Swift (1977). The radiometer output, T H , is plotted on the vertical axis versus the seawater 

temperature, T, on the abscissa with the seawater salinity, S, as a parameter. The Klein and Swift 
model function has been the common choice for e(T,S) in recent remote sensing studies (Le 
Vineet al 1998, Goodberlet et al 1997, Lagerloef 1995). In the remote sensing paradigm, we 
measure and obtain T from an ancillary source. The salinity is then found from Figure 1, at 
the intersection of T B and the seawater temperature T. Given the correct values of T B and T, the 
accuracy of the inversion is limited by the accuracy of the model function, e(T,S). 
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Figure 1 Brightness Temperature vs. Temperature with Salinity as a parameter 


To understand the sensitivity of the seawater salinity inversion to errors in the model function, 
we examine the inversion in more detail. Let us assume that the exact model function is denoted 
by e(T,S) and the model function constructed from measurements is given by £ in (T,S). The error 
in the model function is: Se - e-e . 

nn 

In performing the inversion, it is assumed that the brightness temperature, 7«, and the sea surface 
temperature, T, are known exactly. The model function is then used to find the salinity, S. To 
examine the sensitivity of the algorithm to changes in the model function, we determine the 
change in salinity, SS, for a given change in the model function, Se . It should be noted that the 

dielectric constant is complex and thus & is complex. The results of this calculation are shown 
in Figure 2. Here we have assumed that the measured model function is the Klein and Swift 

model; we have also assumed that the magnitude of e ||e| = V^ /2 +e' r2 j is large (a good 
approximation for seawater at 1.413 GHz). 

Figure 2(a) shows the impact of a 0.2 change in the real and the imaginary parts of the model 
function on the change in the salinity. It illustrates that the error in retrieved salinity can be 
larger than 0.2 psu (the goal for the Aquarius mission) and is larger in cold water than in warm 
water. This should not be surprising since an examination of Figure 1 shows that the constant 
salinity curves come closer together (decreased radiometric sensitivity) as the temperature 
decreases. This fact has been pointed out by Dinnat et al (2003). Figures 2(b) and 2(c) show the 
same curves as in Figure 2(a) except now the change in dielectric constant is completely in the 
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(C) 



Figure 2. Change in Salinity versus Salinity for several temperatures. In (a), (b) and (c) all 
errors have the same magnitude but different real and imaginary parts. 
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real part (b) or in the imaginary part (c). The interesting feature is the effect of changing only the 
imaginary part is much smaller than changing the real part. This says that the real part must be 
measured more accurately than the imaginary part. The details of the sensitivity calculation, used 
to produce these curves, are given in Appendix A. 

2.2 Model Function History 

The model function used in Figure 1 was developed more than 25 years ago by Klein and Swift 
(1977). They used measurements in the frequency range 1-3 GHz that were made by Ho and 
Hall (1973) and Ho et al (1974). These measurements were used to refine the existing 
parameterization for the dependence of £ on T and S at L- band. This parameterization was based 
on the Debye model for the dielectric constant proposed by Stogryn (1971). The result was a set 
of regression curves for these coefficients as a function of T and S that have served the remote 
sensing community well since 1977 and have been used in much of the research demonstrating 
the potential for passive remote sensing of salinity (Blume et al 1978, Goodberlet et al 1997, Le 
Vine et al 1998). Much of this work, however, has been on salinity fields with large gradients, 
such as those occurring in the coastal environment where measurement accuracy of 1 psu is 
adequate. 

The goal of modern remote sensing of the global salinity field from space is almost an order of 
magnitude more ambitious. For example, the goal of Aquarius is 0.2 psu (Koblinsky et al 2003) 
and the goal of SMOS is 0.1 psu (Sivestrin et al 2001). Achieving such accuracy requires careful 
control of all elements of the remote sensing system. This control would include, for example, a 
well designed radiometer that gives accurate values of T«, gives good sources of ancillary data 
such as Sea Surface Temperature and also gives an accurate accounting of environmental factors 
such as emission and attenuation along the propagation path. Modern measurements (Stogryn et 
al 1995, Ellison et al 1998, Meissner and Wentz 2004, Blanch and Agusca 2004; Lang et al 
2008 a, b) suggest that refinements are also needed in the model function. They point out 
that certain aspects of the Ho data may be in error and, as a result, the Klein and Swift model 
function may need modification. The recent measurements of Blanch and Aguasca (2004) 
indicate differences when compared to the Klein and Swift (1977) model function. If valid, these 
differences will result in errors larger than the 0.1 -0.2 psu goal of Aquarius and SMOS. Stogryn 
et al (1995) has suggested an alternate parameterization of the dependence of conductivity on 
salinity and temperature. At L-band, the influence of salinity is primarily through its effect on 
conductivity. The modern definition of salinity, in fact, is based on the measurement of 
conductivity (Lewis, 1980). 

Modern measurements have also been made at higher frequencies (Stogryn 1995, Ellison et al. 
1998, Meissner and Wentz 2004). These measurements also suggest that changes in the model 
function may be needed. Figure 3 shows a comparison of the Ellison et al model function with 
the Klein and Swift model function. The two, based on data in different frequency ranges, 
employ regression methods to fit the salinity and temperature dependence. At first glance, the 
agreement is good (because both employ the Debye model). The scale in Figure 3, however, is 
very large and the differences that appear modest in Figure 3 are far too large to achieve an 
accuracy of 0.2 psu. 

Another example is shown in Figure 12 of section 4.5.2, where measurements made at GW are 
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Salinity = 25 T emperature = 25 C 




Figure 3. Seawater Permittivity from Ellison and Klein and Swift 

plotted together witn tne Klein ana swirt model tor ^ = 3U psu. Aitnougn tne trend with 
temperature is in agreement with the Klein and Swift model, there is clearly an offset (bias), 
especially apparent in the real part of the dielectric constant. The bias is significant for a remote 
sensing goal of 0.2 psu. Similar behavior is seen when comparing the Blanch and Aguasca 
measurements with the Klein and Swift model function. 

The goal of the research work reported here is to make accurate seawater measurements for a 
range of temperatures and salinities characteristic of the open ocean at a frequency of 1.413 GHz 
that is in the center of the Aquarius frequency band. Each measurement at a specific temperature 
and salinity is repeated several times; the mean and variance are also reported. Before describing 
the measurements in section 4, additional details of the background and measurement technique 
are discussed in Section 3. Finally, our conclusions are given in section 5 along with suggestions 
for future work. 

3 Measurement of the Dielectric Constant of Seawater 

In this section, various methods for making measurements of the dielectric constant of seawater 
are described. These descriptions are then followed by a discussion of the results of past 
measurements of seawater. In Section 3.1, a review of the general approaches to the 
measurement process is given along with a discussion of main measurement devices presently in 
use. In Section 3.2, the measurement of the seawater dielectric constant by non GW 
investigators in the frequency band 1400 to 1430 MHz is discussed. Measurements of the 
dielectric constant of seawater at other microwave frequencies will also be reviewed. 
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3.1 General Approach 

The need to know the dielectric constant of solids and liquids at microwave frequencies 
developed with the growth of radar during World War II. The results of this rapid progress have 
been documented by Montgomery (1947) in his chapter on dielectric measurements. Six generic 
measurement methods have evolved: 

• free space reflection and transmission measurements 

• waveguide transmission and reflection measurements 

• reflection and transmission resonant cavity measurements 

These measurement methods are still in use today. Improvements in microwave instruments and 
components however, make today’s measurement results much more accurate . 

The dielectric constant of liquids has been measured by all six methods, although the need to 
confine liquids puts some constraints on the physical configuration. The first method is the free 
space reflection technique that has been employed by Saxton and Lane (1952) and more recently 
by Blue (1980). The method has found use at mm wavelengths. Buckmaster (1990) modified the 
waveguide transmission approach of Nicolson and Ross (1970) to make dielectric constant 
measurements of high loss liquids below 10 GHz. Jenkins et al (1990) have employed a TEM 
coaxial transmission line shorted at one end to make measurements of the dielectric constant of 
several alcohols by the reflection technique. This is a modern version of a reflection technique 
developed by Roberts and von Hippel (1946). Gregory and Clark (2001, 2009) have used 
coaxial transmission and reflection cells to make very accurate measurements of alcohols from 
0.1 to 5.0 GHz. 

In all of the methods except the resonant cavity approach, the frequency can be swept to cover a 
range of frequencies. The resonant cavity method, on the other hand, is usually a single 
frequency measurement. This method, however, is considered the most accurate of the 
techniques (Baker-Jarvis, et al 1998, Baker-Jarvis, et al 2001). There are two types of cavity 
resonators: reflection and transmission. The reflection cavity only has one port and S n is 
measured; the transmission cavity has two ports and S,, is measured. The reflection cavity 

approach is the method that has been used by Ho and Hall (1973) and Ho et al (1974) to measure 
the dielectric constant of seawater. Dielectric resonators can also be used to measure seawater; 
however, their Q may not be high enough (Krupka, 1988), 

Another important advantage of the resonant cavity approach, when making highly accurate 
measurements, is that the seawater can be isolated within the cavity via a glass tube. In 
waveguide transmission and reflection cells, the seawater is in direct contact with the metal walls 
of the cell. Because of the corrosive character of seawater, its interaction with the metal walls 
may lead to serious problems. Stogryn (1995) reported difficulties with repeatability due to this 
problem. To minimize the effects, Stogryn gold-plated the waveguide and washed the walls with 
hydrochloric acid after each measurement. He also mentioned that if the seawater were to 
remain in the transmission cell too long, its salinity would change. 
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3.2 Past Measurements of the Salinity of Seawater 

Over the years, there have been relatively few microwave and mm wave measurements of the 
seawater dielectric constant. Saxton and Lane (1952) conducted one of the earliest investigations. 
They measured the dielectric constant of an NaCl solution that had properties close to those of 
seawater (Ho et al 1974) for frequencies between 9.35 and 48.2 GHz by the surface reflection 
technique. The first measurements of actual seawater were made by Ho and Hall (1973) using a 
reflection cavity technique at 2.65 GHz. This was followed by Ho et al (1974) who reported on 
seawater measurements at 1.43 GHz - again using a reflection cavity technique. Blue (1980) 
made measurements at 97, 103.8, 135.6 and 183.3 GHz using the free space reflection method. 
Stogryn (1995) made measurements of seawater from 7 GHz to 12 GHz by using a waveguide 
transmission cell. Ellison (1998) made measurements from 30 to 105 GHz employing a coaxial 
transmission cell at the lower frequencies and a free space transmission cell at the higher 
frequencies. In 2004 Blanch and Aguasca measured the seawater dielectric constant over a range 
of L band frequencies using a microline transmission cell. In the same year Meissner and Wentz 
(2004) used SSMI satellite measurements to determine, indirectly, the sea surface dielectric 
constant at the SSMI frequencies. Most recently, Lang et al (2008a, b) reported transmission 
cavity measurements made at 1.413 GHz. 

For the measurement of ocean salinity, we are concerned mainly with measurements made in the 
window of 1 .400 to 1 .430 GHz. These include the measurements of Ho et al (1974), Blanch and 
Aguasca (2004) and Lang et al (2008a, b). The other works, however, contain a great deal of 
useful information on measurement techniques. Stogryn et al (1995) in particular conducted very 
accurate measurements by driving down the fluctuations with repeated measurements. 

3.2.1 Measurements made by Ho et al (1974) 

Ho et al (1974) employed the resonant cavity approach to measure the dielectric constant of 
seawater at 1.43GHz. They used the resonant cavity approach because of its known accuracy and 
because the permittivity was needed at only one frequency. They built a cylindrical cavity out of 
copper that was then excited in the TMoio mode. The liquid was then introduced into the cavity 
through a coaxial glass tube with a very small inner diameter. The resulting microwave signal 
was coupled into the cavity by a single loop (reflection cavity approach). The reflection 
coefficient was then measured at the resonant frequency and at some off-resonance frequency 
such as the frequency of the half power point. 

Under the assumption that the amount of seawater introduced into the cavity is small, a 
perturbation theory (Chen et al 2004) can be applied to obtain the following formulas: 

£'-1 = 2 C^-, £' = CA(l/<2) (2a) 

where 

4 /=/,-/. a(i / e)=7;-7r (»> 
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Here £ and are the real and the imaginary parts of the dielectric constant of the sample 
respectively. The quantities f 0 and <2 0 are, respectively, the resonant frequency of the cavity and 

its Q-factor when there is no sample in the cavity. The same quantities without the subscripts are 
the values of the resonant frequency and Q with a sample in the cavity. A derivation of these 
formulas is presented in Appendix B that accounts for cavity wall loss and gives a means of 
computing the error made when using eq(2). 

The quantity C appearing in eq(2a) is called the calibration constant. It can be evaluated by two 
methods: relative and absolute (Ellison 1996). When using the relative technique, another liquid 
of known dielectric constant is measured in the cavity. Since e' and e" are known and 
A/ and A(l/g) are measured, the constant C can be determined from eq(2a). The quantity C can 

also be obtained from the expressions for the electromagnetic fields in the cavity. This is known 
as the absolute calibration technique. It becomes increasingly more accurate as the cavity is 
modeled in more detail. From perturbation theory (Chen et al 2004) it is shown that the 
calibration coefficient appearing in eq(2a) is inversely proportional to the volume of the tube 
holding the sample. For example, a 1% error in tube volume results in a measured dielectric 
constant having a one percent error. 

Ho et al (1974) used both the relative and absolute methods in making their measurements of 
seawater. They employed the relative technique to make the actual seawater measurements. They 
used distilled water as a reference. At that time, however, there was relatively little distilled 
water dielectric data available near their measurement frequency. As a result, they used their 
cavity to measure the dielectric constant of distilled water (absolute method) and used 
expressions for the electromagnetic fields in the cavity to obtain an expression for the calibration 
coefficient. They then carefully measured the volume of the sample in the cavity by two 
techniques. Their findings for the dielectric constant of distilled water have not been checked by 
using another measurement approach at 1 .43 GHz. 

Ho et al (1974) claimed an accuracy of 0.2 % for their measurement of £' and 0.4 % for their 
measurement of ^''for seawater at 1.43 GHz. This is an average accuracy for values obtained 
for temperatures from 5°C to 30°C and for salinities from 5 %o to 35%<f\ But these results 
assume that the dielectric constant of the reference liquid (distilled water) was known exactly. 
However, we have no way of knowing the accuracy of their measurement of distilled water by 
the absolute method (modeling the cavity accurately). They tried to compare their dielectric 
measurements of distilled water with available data but the available data was very sparse. In 
addition, their seawater measurements do not reduce to the distilled water dielectric values when 
these are extrapolated to zero salinity. Ho et al 1974 claimed that this was the effect of organic 
material in the seawater. Ellison (1996), however, has shown experimental results indicating that 
organic material in seawater has little effect on the dielectric constant. Another potential source 
of error is microwave heating of the seawater in the cavity. Swift (2006) mentioned that heating 
effects may have introduced some additional errors in the imaginary part of the dielectric 


** Salinity measurements before 1978 will be reported in parts per thousand (%cj while the measurements after 1978 
will be given in psu. As pointed out in the introduction, the numbers can be used interchangeably in the range of 
salinities and temperatures we are considering. 


9 


constant when the salinity and temperature are high. High salinity values (30 %c to 40 %o) are of 
the most interest in open ocean salinity measurements. 

More recently Stogryn et al. (1995) made measurements of seawater from 7 to 14 GHz by using 
a temperature controlled transmission waveguide of variable length. They measured seawater 
with a salinity range from 0 to 38 psu with temperatures from 0°C to 30°C. They proposed a 
Debye permittivity model with a double relaxation time to better fit the data at higher 
frequencies. In the process of determining the Debye parameters they combined their data (7 to 
14 GHz) with the Ho and Hall (1973) data at 2.65 GHz and with the Ho et al (1974) data at 1 .43 
GHz. They found the deviation from the Debye model for this data set was several percent. 
When they eliminated the Ho data, the error reduced to approximately 0.8%. They conjectured 
that the Ho data contained systematic errors and speculated that these errors could be due to 
microwave heating of the sample within the quartz tube. They pointed out that this was a 
possibility because of the extremely small diameter of the tube. 

In summary, although the Ho et al (1974) measurements were carried out very carefully, their 
accuracy is uncertain, and since the Klein and Swift model function is based in part on Ho‘s 
measurements, the Klein and Swift model function could be in error as well. The need for 
measurements with well documented accuracy to support modern remote sensing - as those 
measurements being performed by SMOS and about to be performed by Aquarius - has led a 
need for new measurements. These include the measurements in support of SMOS by Blanch 
and Aguasca (2004) described below. 

3.2.2. The measurements of Blanch and Aguasca (2004) 

Blanch and Aguasca (2004) have made measurements of the dielectric constant of seawater 
using a transmission cell. They constructed a transmission cell with a stripline structure through 
which seawater flows. Thermometers at the beginning and the end of the cell sensed the water 
temperature in the cell. A large quantity of seawater initially at 0° C was gradually heated as it 
flowed through the transmission cell. During this process, S-parameter measurements were 
continuously taken. The results obtained are close to the Klein and Swift model function for the 
real part of the permittivity. They are somewhat greater than the Klein and Swift model for the 
imaginary part by 1% to 2%. The authors did not give many experimental details such as the 
accuracy of the temperature measurements and error bounds for their measured values. The goal 
of the GW measurements was to provide data with the accuracy needed to achieve the remote 
sensing goals of Aquarius and equally important to document the accuracy with error bars. 

4 Measurements at The George Washington University 

The measurement program at GW that has been supported by this grant began about 4 years ago 
with the goal of focusing specifically on the dielectric constant of seawater at frequencies 
relevant to the remote sensing of salinity,i.e., frequencies in the radio astronomy window at L- 
band centered at 1.413 GHz). A microwave cavity measurement system at 1.413 GHz was 
constructed to measure the dielectric constant of seawater over values of temperature and salinity 
characteristic of the open ocean. The emphasis focused on making the measurements as accurate 
as possible and on documenting error bars for the measurements and calibration. 
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4.1 Measurement System Construction 

After detailed consideration of many techniques available, it was decided that a resonant 
microwave cavity would yield the most accurate measurements given the fact that dielectric 
measurements were needed for only one narrow frequency band (1400 to 1427 MHz). Although 
measurements of alcohols by Gregory and Clark (2009) were made with high accuracy, the 
contact of the salt water with the metal walls for transmission and reflection cells could lead to 
problems. As pointed out previously, Stogryn et al (1995) reported difficulties in repeatability 
due to these problems. The use of a resonant cavity is similar to the method employed by Ho et 
al (1974) with one exception. Ho et al (1974) chose a reflection cavity having only one coupling 
loop to measure Sn- We have chosen a cavity of the transmission type that has two coupling 
loops and measures S 21 . The reason for choosing a cavity of the transmission type is that 
muchless power is coupled into the cavity (Lang et al 2003). The heating problems that may 
have been observed in Ho’s measurements will therefore not be a problem. In addition, 
transmission measurements are inherently more accurate than those of the reflection type. Ho et 
al (1974) used a cavity that resonated at 1430 MHz. This frequency is slightly outside of the 
quiet band. We have chosen a resonance frequency of 1413 MHz that is in the center of the quiet 
band. If one uses the Klein and Swift model function to calculate dielectric constants, the 
difference in dielectric values for the two frequencies is less than 0. 1 %. 

A simplified drawing of the transmission resonant cavity setup is shown in Figure 4. In the 
figure, a network analyzer (NA) is connected to a resonant cavity that is configured to measure 
S 21 over a band of frequencies about the cavity resonant frequency, f„. The subscript, zero, 
indicates that the small glass sample tube shown in Figure 4 is empty. The blue curve generated 




Figure 4. Resonant Microwave Cavity with Network Analyzer 
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on the network analyzer screen and shown to the right of the NA in Figure 4 has a 3db 
bandwidth denoted by BW 0 . The Q of the cavity can be obtained from the 3db bandwidth by 
using the formula: Q„ = fj BW 0 . Seawater is then put into the glass sample tube and the red 
curve is generated. Here the new resonant frequency is /and the Q = f/BW. These quantities can 
be used in eqs(2a) and (2b) to determine the real and imaginary parts of the dielectric constant of 
the liquid in the sample tube. A small reservoir above the sample tube is used to bring the sample 
to the correct temperature before entering the tube. 

The complete experimental setup is shown in Figure 5. In the figure, the resonant cavity is 
shown immersed in a liquid whose temperature is kept constant by the circulator shown on the 
left. The glass superstructure over the cavity is used as a conduit for nitrogen (N 2 ) to apply 
pressure on the seawater to force it into the sample tube. Also shown in this figure are the NA 
and the frequency stabilizer (counter). Detailed photographs of selected parts of the 
experimental apparatus are shown in Figure 6: the cavity with its glass superstructure appears in 
Figure 6a, the inside view of the top cover of the cavity with its microwave coupling loops are 
shown in Figure 6b and the glass sample tube that is placed along the axis of the cylindrical 
cavity and glass reservoir used to hold the seawater before entering the cavity are shown in 
Figure 6c. 

The cavity that has been machined out of naval brass, has a Q of approximately 8000. Brass has 
been used instead of copper since it can be machined more easily to greater precision. The 
cavity can be plated if higher conductivity walls are needed, but that has not been the case up to 
now. The seawater has been introduced into the cavity via a fused quartz tube of outer diameter 
of 3mm. The inner diameter (ID) of the sample tube was chosen to be 0.1mm. This is very small 
but is required if a reasonably large Q (about 4,000) is to be maintained when the seawater enters 
the cavity. This diameter also insures that the accuracy of the perturbation formulas given in 
equations (2) are within 0.1% as pointed out by Ho et al [1974] and demonstrated in Appendix 
B. Because the quartz tube ID is so small, the seawater sample will not usually enter the tube 
unless pressure is applied. To provide a means for the sample to enter the tube under pressure, a 
socket has been attached to the tube (see Figure 6c). Firing the quartz socket to quartz tube has 
proven to be a difficult task since heating the tube tends to close the small hole. The problem has 
been solved by forcing oxygen through the tube under pressure while the tube is fired to the 
socket. The top of the socket is attached to a reservoir to hold the sample. The reservoir has a 
small glass tube within it to hold a thermometer. This thermometer is used to measure the 
temperature of the sample. The reservoir is connected to a supply of dry industrial N 2 . Typically 
2 atmospheres of N 2 were used to push the sample through the tube once it reached the 
temperature of the surrounding liquid. The stable clock is used to keep the NA frequency correct 
to 1Hz. 

The cavity operates in the TM 010 mode. This means that the maximum electric field occurs 
along the axis of the cavity. A potential problem results when the glass tube exits the cavity at 
the top and the bottom. The glass has a low dielectric constant of approximately 4, however, the 
dielectric constant of the seawater sample is many times higher. Ho and Hall (1973) were 
concerned about this problem. Using a perturbation calculation derived by Estin and Bussey 
(1960) they showed that the error introduced was less than 0.1% for a cavity at 2.65 GHz. 
Similar calculations at 1 .4 GHz again showed the effect caused an error of less than 0. 1 %. 
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Figure 5. Experimental setup 

Another consideration is the electrical conductivity between the cavity covers and the cylinder. 
This is a concern since the TMoio mode has currents that go up and down the walls and onto the 
covers. To make a proper connection, silver paste was used along the rim where the walls meet 
the covers. In addition, each cover was screwed down with 16 screws that were torqued to 13 
in-lbs. Tests showed that increasing the pressure did not change the resonant frequency of the 
cavity.The cavity was kept at a stable temperature by immersing it in a bath of water and 
ethylene glycol. A tank was constructed of Vi inch Lexan to hold the cavity and water/ethylene 
glycol 

4.2 Measurement Procedures and Data Acquisition 

The basic microwave measurement device was an HP network analyzer (Model 8722D) operated 
in the transmission mode (measurement of S 21 ). A pair of isolators (17dB backward isolation 
and 0.5 db forward attenuation), one at the input to the cavity and another at the output of the 
cavity were used to reduce the pulling effect of the external circuit on the cavity resonant 
frequency. For some measurements, 5 dB attenuators were used in front of the isolators to 
reduce the return loss of the isolators. For a few measurements, a 17 dB RF low noise amplifier 
(LNA) was used to increase the signal to noise ratio of the output signal. The amplifier was later 
eliminated since gain variations due to room temperature changes caused S 21 to vary. This setup 
is shown in Figure 7. In Appendix C.l, tables show which measurement configurations were 
used for different measurements. 
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Figure 6 (a) -cavity with glass reservoir superstructure; (b) - cover of resonator showing 
location of coupling loops; (c) left side - thin quartz tube that holds seawater inside cavity 
together with quartz coupling socket; (c) right side - reservoir with interior glass tube to hold 
thermometer. 
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The microwave power is coupled in and out of the cavity via coupling loops. These loops can be 
raised or lowered to decrease or increase the coupling into the cavity. Normally, a cavity in the 
transmission mode is operated with very weak coupling. Under these conditions, the unloaded Q 
of the cavity is approximately equal to the loaded Q. It is the loaded Q that is measured and the 
unloaded Q that is used in eq (2a) to find the imaginary part of the dielectric constant. This 
topic is discussed in more detail in Appendix D. There, it is shown that the unloaded Q does not 
necessarily have to be very close to the loaded Q in order to be used in eq(2a). In the appendix, 
the A(l/0 used in eq(2a) is calculated for the unloaded and loaded cases. For all measurements 
made in this report, the two are virtually the same even though in some cases, the loaded and 
unloaded Q differ substantially. 

The NA was set to a center frequency of 1413 MHz with a span of 50KHz. As measurements 
spanned the temperature range from 10°C to 35°C, the resonant frequency varied about 700 KHz. 
As a result, as the temperature changed, the center frequency of the sweep had to be adjusted 
slightly to keep the resonant curve on the NA screen. The analyzer was set for 8 averages, 1601 
points, an IF bandwidth of 100Hz and a power output of -5 dBm. This is the maximum output 
for this NA model. At the beginning of each measurement, a thru-calibration was done between 
points a and b as shown in Figure 7. 

Temperature was measured using three thermistors purchased from US Sensor with an accuracy 
of 0.005°C. A rectangular block of naval brass was mounted to the top and bottom of the cavity. 
A thermistor was tightly fitted into a cylindrical hole in each block. The third thermistor was 
inserted in the small glass tube inside the reservoir to measure the temperature of the seawater 
sample. In Figure 6(a), the thermometer is shown inserted into the reservoir tube. 

Measurements were made to determine how long it would take distilled water, at room 
temperature, to reach the temperature of the surrounding fluid after it was introduced into the 
reservoir. To do this, the third thermistor was placed inside the reservoir tube with some alcohol 
so that the thermistor made thermal contact with the glass. Because of the length of the 
thermistor (3 inches), a large amount of distilled water had to be used to completely cover the 
length of the thermistor. The measurements showed that, for all temperatures in the range 
(10°C<T<35°C) 20 minutes was more then enough time for the distilled water to come to 
equilibrium with the surrounding liquid, These measurements were taken using distilled water 
rather than seawater because calibrated seawater is quite expensive (approximately $150 for a 
small ampoule) and the amount of seawater needed was quite large. After the initial tests with 
distilled water, the third thermistor was used to measure the room temperature. Inside the fixed 
temperature bath, repeated measurements showed that the temperature difference between the 
thermister on top of the cavity and the one on the bottom differed by less than 0.005° C. Toward 
the end of the measurement period, the lead to the top-mounted thermister broke off flush with 
the thermister head making it difficult to repair. The measurements proceeded with only two 
thermisters: one mounted on the bottom of the cavity and one placed to measure the room 
temperature. 

After the NA thru-calibration is completed, the Visual Basic program is started. Its measurement 
sequence is as follows: the NA performs 8 sweeps across 1601 frequency points. An 8 sweep 
average of S 21 is computed for each frequency point. By using the internal functions of the NA. 
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Figure 7. Measurement Setup 

the center frequency, the 3db bandwidth and the Q of the resonant curve are computed. These 
values are only computed if the resonant curve, including the 3db points, are on the NA screen. 
Once these calculations are complete, the Visual Basic program downloads the data into a single 
column of an Excel spread sheet. First, the start and stop time of the measurements are 
downloaded, then the three temperature measurements at the start time and the three 
temperatures at the stop time are downloaded. This is followed by the center frequency, 3dB 
bandwidth, the Q and the magnitude of S 21 in dB at the resonance frequency. Following these 
quantities, the 1601 values of the magnitude of S 21 in dB at each frequency are downloaded. 
These 1601 points can be used to independently determine the resonant frequency and bandwidth 
of the resonant curve, if desired. 

The generation of one column of data represents one measurement cycle. The cycle takes a little 
less than 3 minutes to complete. Once it is completed, the Visual Basic program starts the next 
cycle. It is downloaded into the second column of the spread sheet. The process keeps on going 
until it is stopped. 

A typical seawater experiment consists of many cycles. First, the cavity is brought to the desired 
temperature and allowed to stabilize. Next, 5 to 10 empty measurements, i.e., no seawater in the 
tube, are made. The seawater is then put into the glassware by opening the glass cap at the top of 
the glass structure (see Figure 6(a)) and inserting a small nylon tube down to the top of the 
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cavity. The nylon tube is used so that the walls of the glassware won’t be contaminated with 
seawater. Only the tube and the very bottom of the glassware need to be washed after the 
experiment is finished. The seawater is inserted into the tube by using a syringe with a 50 jum 
filter on it. The glass cap is replaced, and the seawater is allowed to stand about 20 minutes until 
it reaches the temperature of the surrounding liquid. Because of the viscous nature of seawater, 
it usually does not enter the tube. To start the measurements, N 2 ( at a pressure of about 2 
atmospheres) is used to push the seawater into the sample tube. The pressure is kept constant 
until all the seawater leaves the reservoir and the tube is empty. Next the pressure is released 
and distilled water is added to clean the cavity. Once the distilled water is added, the pressure is 
turned back on. Enough distilled water is added so as to clean the glass surface that has been in 
contact with the seawater. This cleaning procedure can take from 15 to 30 minutes depending on 
how much distilled water is added. The last step is a set of empty tube measurements to confirm 
that the cavity has indeed returned to its empty state. The whole procedure to measure the 
seawater at one temperature and one salinity can take from 1 to 3 hours. It should be pointed out 
that, after the tube is cleaned with distilled water, the resonant frequency does not usually return 
to exactly the same value it had when the experiment was started a few hours earlier. The 
resonant frequency of the cavity exhibits a slow drift over time. This drift is discussed in 
Appendix C.3. 

An Excel worksheet has been coded to plot measured variables as a function of time. Some data 
from an experiment performed on May 19, 2008 to measure the dielectric constant of seawater 
having a salinity of 34.999 psu and temperature of 25°C is shown in Figures 8 and 9. In Figure 
8, the resonant frequency and the Q are shown as a function of time. Each diamond on the plot 
represents the value of the plotted variable in a particular column of the spread sheet. The time 
scale is expressed in terms of the number of measurement cycles. Thus if each experiment takes 
approximately 3 minutes, the value 10 on the time scale represents 30 minutes. The total time 
shown for the experiment represents about 3Vz hours (3 x 70 = 210 minutes). In Figure 9, the 
cavity temperature and the room temperature are shown. We see that over the total time required 
for the measurements (3*/2 hours), the cavity temperature remained stable at 25.04+0.01 5°C while 
the room temperature varied. The wave like character of the cavity temperature graph reflects 
the cyclic switching of the circulator between cooling and heating. 

Consider the measurements of Figure 8 again. At the beginning, the sample tube in the cavity is 
empty. An atmosphere or two of N 2 is applied to ensure that the sample tube is empty at the 
start of the experiment. As a result, the Q remains about 7800 for the first 18 cycles. At 
approximately the 10 lh or 1 l lh cycle, the seawater is added to the reservoir, however it does not 
enter the tube immediately since it is too viscous. After waiting 7 to 8 cycles for the seawater to 
reach 25°C (about 20 minutes), a pressure of 2 atmospheres of N 2 is applied. It pushes the 
seawater into the tube and the Q falls to about 4600. Now 6 measurement cycles occur before all 
the seawater is pushed through the tube. At this time, the pressure is released and distilled water 
is added. The Q levels out at slightly less than the free space Q since distilled water has a very 
low loss. After about 12 cycles the distilled water is completely expelled from the reservoir and 
the tube. The Q returns to the empty state. If it does not, a second washing with distilled water 
is made. 
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Figure 8. Measurements of the resonant frequency and Q as a function of time. 
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Figure 9. Measurements of cavity and room temperature as a function of time. 
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To compute the real part of the dielectric constant of seawater, we will need the change in the 
resonant frequency when the tube is empty and when it has seawater in it. The change in 1/Q 
will be needed to compute the imaginary part of the dielectric constant. Three points are selected 
from the empty measurements right before the transition to the seawater state. These are 
recorded as sample 1, 2 and 3 in Table 1 under Empty Measurements. Three points 
(measurement cycles) are selected from the 6 cycles when the seawater is in the tube and entered 
into Table 1 under Seawater Measurements. The cycles closest to the transition point are usually 
chosen. Finally three points from the empty portion - after the distilled water has been expelled 
from the tube - are chosen and entered into the table under Ending Empty Measurements. In 
all three cases, the three points chosen are contiguous points. The average of the three samples 
is computed in each case. These averages are the values that will be used in the calculation. The 
standard deviation is also found. This provides an important check. If it is too large, other points 
will be chosen and if this is not satisfactory the measurement will be discarded. 


EXPERIMENTAL DATA 


EMPTY MEASUREMENTS 


May 19,08 

Res. Freq (Hz) 

Q 

Cavity T 

Room T 

Time 

Sample 1 

1,413,076,041 

7,827 

25.05 

25.68 

11:47:10 

Sample 2 

1,413,075,866 

7,843 

25.05 

25.80 

1 1 :49:59 

Sample 3 

1,413,075,964 

7,832 

25.06 

25.88 

1 1 :52:47 

Average 

1,413,075,957 

7,834 

25.06 

25.79 

1 1 :49:59 

Stdev 

88 

8 

0.00 

0.10 


SEAWATER 

MEASUREMENTS (34.997psu) 



May 19,08 

Res. Freq (Hz) 

Q 

Cavity T 

Room T 

Time 

Sample 1 

1,413,015,075 

4,581 

25.05 

26.03 

12:01:14 

Sample 2 

1,413,014,899 

4,593 

25.05 

26.06 

12:04:03 

Sample 3 

1,413,014,980 

4,599 

25.06 

26.09 

12:06:51 

Average 

1,413,014,985 

4,591 

25.05 

26.06 

12:04:03 

Stdev 

88 

9 

0.00 

0.03 


ENDING EMPTY MEASUREMENTS ( 

Once Distilled Water Cleaning) 

May 19,08 

Res. Freq (Hz) 

Q 

Cavity T 

Room T 

Time 

Sample 1 

1,413,076,114 

7,833 

25.05 

26.09 

12:54:38 

Sample 2 

1,413,076,208 

7,834 

25.05 

26.22 

12:57:27 

Sample 3 

1,413,076,020 

7,840 

25.05 

26.36 

13:00:16 

Average 

1,413,076,114 

7,835 

25.05 

26.23 

12:57:27 

Stdev 

93 

4 

0.00 

0.14 



Table 1. Experiment data from a measurement made on May 19, 2008 of 
seawater having salinity of 34.999 psu at a temperature of 25°C. 
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4.3 Validity of Perturbation Theory 

The perturbation equations presented in eqs(2a) and (2b) have existed in the literature since 
resonant cavity measurements were first used to measure dielectric constants more than fifty 
years ago. In Chen et al (2004), the usual derivation of the formulas is given. This derivation 
has several deficiencies. First, it assumes that the empty cavity loss can be neglected, and 
second, it doesn’t indicate how to compute an error term. In appendix B, a derivation of the 
perturbation equations is presented that accounts for the empty cavity losses. It also gives first 
order error terms that can be computed. 

The perturbation equations derived in Appendix B are identical to those given in eqs.(2a) and 
(2b). They are: 


e'-\ = 2C — 

f 


e" = CM\IQ) 


(3) 


The only difference between the above equations and those given in eqs(2a) and (2b) is the 
calibration constant C. In Appendix B, C is given as a product 


C=ZC 0 (4) 

where Co is a constant of the system while X depends on the sample properties and the 
temperature. The constant Co is the ratio of the average energy in the empty cavity to the 
average energy in the empty tube. It is a constant of the system; it only changes slightly as the 
temperature changes and the cavity expands or contracts a slight amount. The multiplicative 
factory, on the other hand, is an energy storage factor which changes with the type of liquid 
being measured and the temperature. If the perturbation is small enough thenZ is close to 1 and 
the calibration coefficient is close to Co, a constant of the system. In Appendix B, an expression 
for x * s given in terms of the fields in the cavity when it is empty and when it has a liquid in it. 
By solving the boundary value problem for the fields in the cavity under both these conditions, 
the energy factor % can be computed. It is seen that as the inner tube radius becomes smaller 
(smaller sample size) the factor X goes to 1 . 

4.4 Calibration 

Before seawater dielectric measurements can be made, the calibration coefficient, C, must be 
known. From eq(3), it is seen that C depends on the dielectric constant of the sample itself 
unless X~\. Therefore, to make the formulas useable, it must be required thatZ~l. As has 
been pointed out previously, X can be made close to 1 by making the tube radius small enough. 
More discussion on this point occurs later in this section. 

As was mentioned earlier in the report, the value of C can be determined by two calibration 
methods: absolute and relative calibration. For both methods the assumption is made that C ~ 
Co. For the absolute calibration method Co is calculated from its definition that is given in 
Appendix B. This requires detailed electromagnetic modeling of the cavity and accurate 
knowledge of the sample tube ID. The relative approach uses another reference liquid whose 
dielectric constant is accurately known a priori. The method is quite accurate since it is less 
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dependent on any imperfections in the cavity and it does not require accurate knowledge of the 
sample tube ID. This method, however, depends on an accurate knowledge of the reference 
liquid’s dielectric constant. Since accurate measurements of methanol’s dielectric constant for 
the temperature range and the approximate frequency of interest are available, the relative 
approach has been used in these measurements. 

4.4.1 Calibration Liquids 

To use the relative calibration technique, the complex dielectric constant of some liquid must be 
accurately known. Consequently, the total accuracy of our measurements will depend on the 
accuracy of the seawater measurements, as well as, on the accuracy of the calibration coefficient. 
The accuracy of the calibration coefficient, in turn, depends on the accurate knowledge of the 
dielectric constant of the reference liquid. As will be demonstrated, it is also necessary to have a 
liquid whose dielectric constant is not too different from the dielectric constant of seawater. 

Two easily obtainable liquids that meet these criteria are methanol and distilled water. Gregory 
and Clark (2009) have made accurate measurements of methanol at 20 individual frequencies 
from 0.1 to 5.0 GHz for a temperature range from 10°C to 50°C in increments of 5 degrees. 
Their measurements have an accuracy of approximately 0.3% in most cases of interest to us. 
The measurements have been made with coaxial transmission lines of different lengths and a 
coaxial reflection device. In all cases each result has been checked using two measurement 
devices. The other liquid is distilled water. Kaatze (1989) has made measurements of distilled 
water over a frequency range from 1.1 GHz to 57 GHz for various temperatures in our range by 
using a microwave double beam interferometer. From these values he has constructed a model 
function. He claims an accuracy of better than 1 % for the real part and 3% for the imaginary part 
in the frequency range from 1.1 GHz to 5.3 GHz. Since the accuracy of the methanol 
measurements is better than the distilled water measurements, we will use methanol as our 
calibration liquid. Distilled water, however, does have an advantage since its dielectric constant 
is closer to the dielectric constant of seawater. 

4.4.2 Variation in Calibration Coefficient with Liquid Type 

In this subsection, we develop the method for the use of a reference liquid to calibrate the sample 
equations and the dependence of this procedure on the choice of the calibration liquid. The 
calibration coefficient, C, appearing in eq(4) is a product of C 0 and X , a factor that depends on 
the dielectric constant of the liquid being measured. In relating the calibration measurements to 
the seawater measurements, there are two systems of equations to consider: the sample equations 
and the reference equations. The sample equations are: 



= c,a 0/2 * > 


(5) 


and the equations for the reference liquid are: 



fr 


£ r = C r A(1 /Q r ). 


(6) 
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Sample tube diameter (mm) 


Figure 10. Plot of Error(%) = 100x(^ -1 ) versus the inner diameter of the sample tube. 


By employing eq(4), the calibration constants are each related to C 0 as follows: 

c, = Z,C 0 C, = Z r C„ ( 7 ) 

Combining these the two equations above gives 

C, = (Z./Z,)C r (8) 

To determine the calibration coefficient, C s , for use in the sample measurement equations, the 
following procedure is followed: first, use either the known values of e r or e r in eq(6) to 

determine C r and second, assuming X S IX, is known, compute C.,from eq(8). The ratio of then’s 
is not known but it can be made close to 1. In practice, this is done by making the inner tube 
diameter small enough and choosing a reference liquid with a dielectric constant close to that of 
the sample. An approximate expression for x is given in Appendix B. Plots of this ratio as a 
function of tube radius with methanol as the reference liquid at 20° C ( e r = 27.96 + il 1.20) and 
with distilled water as a reference liquid at 20° C ( e r = 79.75 + j'5.83) are shown in Figure 10. 
The sample dielectric constant has been assigned a large value (f v = 100 + /1 00) to assume the 
worst possible case. For real seawater samples, the corresponding error will be slightly smaller. 
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An examination of Figure 10 shows that for a tube with an inner diameter of 0.1mm, the error in 
approximating C s by C r for methanol is about 0.06% which is negligible. From this point on, we 
will refer to a single calibration coefficient C, determined from measurements of methanol as 
reference liquid. 

It should be pointed out that if X s I Z r is not chosen close to 1, then the measurement procedure 
will be a function of the dielectric constant being measured making the value of /Z, different 
for each value of £ s - not a good situation. Note from Figure 1 0 that when distilled water is used, 
a larger tube radius can be employed for same error. Normally distilled water would have been 
chosen for a calibration liquid but values of its dielectric constant at the frequency and 
temperatures of interest are not known with the same accuracy as for methanol. 

4.4.3 Calibration Measurements Using Methanol 

Multiple measurements of methanol were made over the period from January 26, 2008 to May 
13, 2008 to determine a calibration coefficient for tube #12. Table C. 1-1 in Appendix C.l shows 
the microwave network (see Figure 7) that was used for each measurement. The table indicates 
that isolators were always employed. For the data taken in April and May of 2008, 5 dB 
attenuators were also used. These microwave devices were used to isolate the cavity from its 
external environment. The strategy was not completely successful in stabilizing the long term 
frequency drift, possibly due to the fact that the temperatures of the attenuators and the isolators 
also changed as the room temperature varied. 

The measurements of resonant frequency and Q for the empty cavity and for the cavity with 
methanol in it are shown in Table 2. The methanol used was purchased in early January, 2008 
from Fisher Scientific. Its product number was A935-4 and it was rated as having low water 
content. In Table 2, the difference frequency, A f r , and the difference A(l/<2 r ) are computed for 
use in eq(6). The calibration coefficient, C can be calculated using either the equation 
determining e' r or the equation determining £* in eq(6). We will use both equations and obtain 

two values of C. They are referred to in the Table 2 as Cl (determined from the £ r equation) and 
C2 (determined from the £* equation). The closeness of Cl and C2 will be used as a measure of 
the self consistency of the measurement. 

Before calculating Cl and C2, £ r and £* must be determined from the tables of Gregory and 

Clarke (2009). Since the tables do not provide the dielectric constant for a frequency of 1.413 
GHz, interpolation will be used. This process is performed by sixth degree polynomial 
interpolation, and the results are given in the caption to Figure 11. Now that the complex 
dielectric constant of methanol at 20° C has been determined, the fonnulas in eq(6) can be used 
to determine Cl and C2. The results of all the methanol measurements that have been made are 
listed in Table 2. From these experimental results, the values of Cl and C2 are also computed 
and listed in the table. 

Rather than use all the data in Table 2 to determine C, the percent difference between Cl and C2 
will be used to eliminate some data points as being too inconsistent. The quantity 100x(Cl/C2-l) 
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Date 

Empty f 0 (Hz) 

Methanol 

fr(Hz) 

Af r (Hz) 

Empty 

Q 0 

Methanol 

Qr 

A(1/Q r ) 

Cl 

C2 

|C1/C2-1 1 
xlOO 

01/26/08 

1,413,279,303 

1,413,255,197 

24,107 

7,741 

6,979 

1 .4094E-05 

790,280 

794,664 

0.552% 

01/26/08 

1,413,278,552 

1 ,413,254,151 

24,402 

7,752 

6,981 

1.4241E-05 

780,725 

786,462 

0.729% 

01/26/08 

1,413,278,164 

1,413,254,037 

24,127 

7,766 

6,993 

1.4219E-05 

789,613 

787,701 

0.243% 

01/27/08 

1,413,275,110 

1,413,251,056 

24,054 

7,779 

7,004 

1.4226E-05 

792,018 

787,270 

0.603% 

01/31/08 

1,413,258,135 

1,413,233,735 

24,400 

7,692 

6,938 

1.4132E-05 

780,757 

792,553 

1 .488% 

01/31/08 

1,413,257,907 

1,413,233,909 

23,998 

7,640 

6,925 

1.3516E-05 

793,846 

828,648 

4.200% 

02/01/08 

1,413,265,777 

1,413,241,410 

24,366 

7,675 

6,922 

1.4180E-05 

781,850 

789,864 

1.015% 

02/01/08 

1,413,268,742 

1,413,243,084 

25,657 

7,658 

6,914 

1.4062E-05 

742,521 

796,473 

6.774% 

02/04/08 

1,413,260,615 

1,413,236,397 

24,218 

7,632 

6,899 

1.3918E-05 

786,636 

804,696 

2.244% 

04/1 0/08 

1,413,218,048 

1,413,193,860 

24,188 

7,375 

6,673 

1.4261 E-05 

787,588 

785,381 

0.281% 

04/19/08 

1,413,208,681 

1,413,185,479 

23,202 

7,883 

7,103 

1.3937E-05 

821,052 

803,628 

2.168% 

04/24/08 

1,413,192,694 

1,413,169,244 

23,451 

7,870 

7,086 

1.4061 E-05 

812,325 

796,529 

1.983% 

04/24/08 

1,413,193,573 

1,413,170,256 

23,317 

7,865 

7,080 

1.4096E-05 

816,994 

794,529 

2.827% 

04/25/08 

1,413,194,693 

1,413,171,086 

23,607 

7,869 

7,080 

1.41 71 E-05 

806,958 

790,369 

2.099% 

04/29/08 

1,413,226,919 

1,413,202,506 

24,414 

7,633 

6,891 

1.4109E-05 

780,302 

793,842 

1.706% 

04/30/08 

1,413,223,373 

1,413,199,184 

24,189 

7,569 

6,838 

1.4120E-05 

787,558 

793,224 

0.714% 

04/30/08 

1,413,223,731 

1,413,199,468 

24,263 

7,573 

6,833 

1 .4299E-05 

785,157 

783,255 

0.243% 

05/1 3/08 

1,413,207,779 

1,413,183,880 

23,899 

7,846 

7,068 

1.4039E-05 

797,106 

797,762 

0.082% 

05/1 3/08 

1,413,205,702 

1,413,181,833 

23,870 

7,852 

7,074 

1.4003E-05 

798,085 

799,832 

0.218% 

Average 

1,413,237,763 

1,413,213,672 

24,091 

7,719 

6,962 

1.4089E-05 

791,125 

795,088 



Table 2. Chart of measured methanol data and calculated values of Cl and C2 for each data 
point. All measured data with 100x(Cl/C2-l) < 1.5% are marked in yellow. 


Calibration Analysis, Cl 

Average 

788,249 

Stdev 

5996 

Error (%) 

0.76% 


Calibration Analysis, C2 

Average 

791208 

Stdev 

5333 

Error (%) 

0.67% 


Table 3. The charts above give the average and standard deviation of Cl (Table 3a) and of C2 
(Table 3b) over measurements that have a consistency of better than 1.5%. The error(%) 
=(Stdev/ Average) xlOO. 

is computed in the last column. Only those experiments having C’s that differ by less than 1 .5% 
will be used in the calculation of the average C and its standard deviation.. These experiments are 
shown in yellow in Table 2. The average values of Cl and C2 over measurement values that 
have a consistency of at least 1 .5% are given in Table 3. Also shown in the tables is the standard 
deviation over the consistent measurement points. The value of <C/> will be used as C i.e., C = 
<C1>. 

4.5 Seawater Dielectric Measurement Data 

In this section, the measurement of seawater samples will be discussed. The section first 
addresses where the seawater samples were obtained. Following this, the measurement 
configuration and times the measurements were taken are discussed. Next, the actual data with 
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Methanol at 20C Dielectric Constant vs Frequency 

35 — 



0.00 0.50 1.00 1.50 2.00 2.50 3.00 3.50 

Frequency (GHz) 

-*-*■ y = -0.1354x + 1.4489x - 6.3724x + 15.296x -21.71 lx + 1 1 ,031x + 30.892 
”*~ I "y=-0 0H9x +0.0758x - 0.0066x‘ - 0.6007x - 1.2626x + 10. 893x - 0.2845 

Figure 1 1 . Polynomial fits to the Gregory and Clarke data for methanol at 20° C. At a 
frequency of 1 .413 GHz the interpolated values are: £ r - 27.96 and e"= 1 1.20. 


its means and variances are presented. This is followed by an analysis of the errors. 

4.5.1 Seawater Samples 

The seawater was purchased from Ocean Scientific International Limited (OSIL) located in 
Havart, Hampshire in the United Kingdom. They sell IAPSO Standard Seawater. The standard 
seawater as approved by the International Association for Physical Sciences of the Ocean 
(IAPSO) is the only transfer standard for Practical Salinity Scale that is recognized by the major 
oceanographic bodies. 

The samples are prepared from surface seawater collected from the North Atlantic and then 
pumped through filters (0.2 mm) into a PVC lined tank and thoroughly mixed. Once purified by 
filtration and UV radiation, distilled water is added to adjust the salinity and the water is sealed 
in pre-washed glass bottles (capacity 200 cm’). The ratio of electrical conductivity of each batch 
of seawater to the conductivity of a defined standard potassium chloride (KC1) solution is used to 
measure the salinity in psu. This is done by using a high precision salinometer. 

The seawater is sold in several standard packages. GW purchased the linearity pack. It consists 
of 10 bottles (200ml each) of seawater having the following salinities: 4 bottles of 35 psu, 2 
bottles of 38 psu, 2 bottles of 30 psu and 2 bottles of 10 psu. Two seawater linearity packs have 
been purchased: one in November of 2002 and another in March, 2008. OSIL says that the 
seawater has negligible change in salinity in a shelf life of 3 years. Bacon et al (2007) have 
shown that the salinity changes by less than .001 in 282 weeks. Only seawater purchased in 
November, 2002 has been used in the experiments reported here. 
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4.5.2 Seawater Measurements 

The seawater dielectric constant measurements reported here have been made over the period 
from December, 2007 to June, 2008. Seawater having salinities of 30.002, 34.997 and 38.274 
psu was measured. The seawater having salinity of 30.002 psu was measured during the period 
from December, 2007 to February, 2008; the 34.997 psu seawater was measured from May to 
June 2008 and the seawater of salinity 38.274 psu was measured from March to May, 2008. 

During this time, only one quartz tube (#12) was employed and thus the methanol calibration 
measurements can be applied to all measurements. As noted previously, calibration 
measurements were made from January, 2008 to May, 2008. As more calibration data became 
available, the calibration coefficient was recomputed. As a result, measurement data appearing in 
Lang et al (2008a, 2008b) may differ slightly from the results given in this report. In addition, 
other seawater measurements have also been added to the data set since the publication of those 
papers. Although the curves look slightly different than the ones appearing in the conference 
papers, the conclusions published are essentially the same. 

The configuration of the microwave network connecting the NA to the cavity has been changed 
during the experiments. The network configuration for each seawater measurement made is 
shown in Appendix C.l in Tables C.l-2 for seawater of salinity of 30.002 psu, in Table C.1-3 for 
seawater of salinity of 34.997 and in Table C. 1 -4 for seawater of salinity of 38.274 psu. 

The measurements of 30 psu seawater started in December of 2007. At that time, two isolators 
were employed at the input and output of the cavity. Previous experimentation had shown that 
the long term stability of the resonant frequency measurements was less influenced by external 
temperature changes when isolators were used. A 17 dB amplifier was also employed to 
increase the signal-to-noise ratio. After the measurements were made in February, 2008 it was 
decided to eliminate the amplifier. Since the amplifier was in the calibration loop, changes in 
amplifier gain caused by room temperature changes, affected the NA calibration. The change in 
calibration did not affect the cavity resonant frequency or the Q. The resulting measurements 
made while the amplifier had been in the circuit were accurate. 

Measurements made during March and April of 2008 employed only the isolators. Before the 
May measurements began attenuators were added in front of the isolators to reduce their return 
loss. The 34.997 psu measurements employed isolators and attenuators since they were made in 
May and June, 2008. The 38.274 psu measurements were made in March and in May, 2008. The 
March measurements used only isolators while the May measurements used isolators and 
attenuators. 

The individual seawater measurements are given in Appendix C.2 in 3 tables labeled Table C.2-1 
for the 30.002 psu data, Table C.2-2 for the 34.997 psu data and Table C.2-3 for the 38.274 data. 
In these tables the date for each measurement is given along with the change in resonant 
frequency, Af, and the change in 1 /Q, A(l/<2). Using these quantities together with the 
calibration constant C (C = 788,249) in eq(2a) and eq(2b), the real and imaginary parts of the 
complex dielectric constant of seawater are calculated. They are compared with the values 
predicted by the Klein and Swift (KS) model function for seawater in the last two columns. It 
should be pointed out that in making the calculation of the imaginary part, the measured Q is 
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used. It is actually the unloaded Q that should have been used. In Appendix D, it is shown for a 
representative subset of measurements that there is virtually no difference between A(l/0 
evaluated with the loaded Q or the unloaded Q. Although the loaded and unloaded Q values 
differ from each other, the difference cancels out when computing A(l/0. 

An examination of the seawater measurement tables in Appendix C.2 shows that the change in 
frequency, A/, lies between 58 to 67 KHz. These frequency differences must be accurate to three 
places to give real dielectric values that have errors of less than 0.2%. The data also shows that 
the empty Q is approximately 7800 while the Q’s with the seawater vary from 4000 to 5000 
depending on the salinity of the seawater. Similarly, the Q must be accurate to three places to 
give the required accuracy. The means and standard deviations for each measurement set at a 
given salinity and temperature are presented in Table 4. The mean values at a given salinity are 
then used to construct graphs of real and imaginary parts of the dielectric constant as a function 
of temperature for a given salinity. The standard deviation of each measurement point is shown 
by a vertical bar that goes up by <7 and down by o where <7 is the standard deviation. The KS 
model function has been drawn on the same plots for comparison purposes. These plots are 
shown in Figure 12 for salinity of 30.002 psu, in Figure 13 for salinity of 34.997 psu, and in 
Figure 14 for salinity of 38.274 psu. It should be pointed out that an average of 3.3 points per 
measurement has been used to compute these statistics. 


GWU Seawater Data 

Temp (Celsius) 

psu 

e' (Mean) 

e' (Standard Deviation) 

e" (Mean) 

e" (Standard Deviation) 

10 

30.002 

75.01 

0.47 

49.66 

0.19 

15 

30.002 

73.08 

0.31 

53.90 

0.55 

20 

30.002 

72.11 

0.49 

58.65 

0.34 

25 

30.002 

70.92 

0.83 

63.79 

0.49 

30 

30.002 

69.32 

0.45 

68.24 

0.51 

35 

30.002 

68.68 

1.45 

73.71 

0.36 

10 

34.997 

74.20 

0.13 

55.06 

0.04 

15 

34.997 

71.98 

0.08 

60,35 

0.14 

20 

34.997 

70.64 

0.00 

65.92 

0.06 

25 

34.997 

69.08 

0.15 

71.33 

0.08 

30 

34.997 

68.44 

1.01 

76.57 

0.90 

35 

34.997 

67.19 

0.20 

82.51 

0.16 

10 

38.274 

73.29 

0.14 

59.26 

0.21 

15 

38.274 

71.36 

0.71 

64.80 

0.53 

20 

38.274 

70.17 

0.42 

69.82 

0.13 

25 

38.274 

69.03 

0.35 

76.46 

0.26 

30 

38.274 

67.62 

0.52 

82.06 

0.86 

35 

38.274 

66.54 

0.97 

88.14 

1.78 


Table 4. Means and standard deviations for complex dielectric constant measurements of 
seawater for each salinity and temperature. 
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An examination of the six graphs shows that the GW data curves fall a little below the KS model 
results. To assess the amount more accurately, the % differences are calculated in Table 5. The 
table values show that on the average the GW data is about 1 .5% below the KS curves for the 
real part of the dielectric constant. The imaginary curves are overall about 2% below the KS 
data. The imaginary data does show an increasing error as the salinity and temperature increases. 
Looking at the individual curves, the KS curves appear very smooth since they are generated 
from algebraic expressions whose coefficients have been estimated from various data sources. 
The GW curves, on the other hand, are not as smooth since straight lines have been used to 
connect the mean measurement values. Regarding the standard deviations, it appears from 
examining all six graphs that the crincreases with increasing temperature. 

4.5.3 Uncertainty in Seawater Measurements 

In this section, the uncertainty in the seawater measurements will be quantified. The uncertainty 
comes from three major sources: 

• uncertainty in the seawater measurements themselves (Repeated measurements do 
not give exactly the same value of the complex dielectric constant each time they are 
performed.) 

• uncertainty in the measurement of the calibration constant (There are 11 sets of 
acceptable measurements of methanol in Table 2. Each gives a slightly different 
value of the calibration constant, C .) 

• uncertainty in the reference complex dielectric constant of methanol (The 
measurements reported by Gregory and Clarke (2009) have an average uncertainty of 
approximately 0.05 for the range of dielectric constants considered in this report.) 

The cumulative effect of all these uncertainties will be computed. Since there are only a few data 
points for each measurement at a particular value of salinity and temperature, all measurements 
at the same temperature will be combined in making these calculations. An uncertainty will be 
found for the real and imaginary parts of the dielectric constants at each of the six temperatures 
measured. 

Consider the uncertainty in the real dielectric constant first. The expression for the real dielectric 
constant from eq(2a) is: 


e' = 2CAf/f 0 (9) 

The value of £'can be considered a random variable that depends on the product of two random 
variables C and A f . The variation of f a is very small and can be neglected. If it is assumed that 
the random variables C and A f are independent then the following expression is obtained for the 
mean of £'\ 

, . 2(C)(A/) 

(g)= J “ (10) 

J o 

The fluctuations for £' can be expressed as Se' = £'-(£'} . Now taking the fluctuations of eq(9) 
gives 
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Figure 12. Real and imaginary dielectric constant of 30.002 psu seawater; GW 
results compared with KS model function. 
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Figure 13. Real and imaginary dielectric constant of 34.999 psu seawater; GW 
results compared with KS model function. 
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Dielectric Constant (Real) vs. Temperature for Seawater at 38.274psu 



— GWU 

Klein-Swift 



Figure 14. Real and imaginary dielectric constant of 38.274 psu seawater; GW 
results compared with KS model function. 
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Mean Measured Dielectric Constants vs. KS 


psu 

Temp°C 

<£’ > 

e’(K&S) 

£’%diff 

A 

V 

£"(K&S) 

£"%diff 


30.002 

10 

75.01 

75.97 

1.28 

49.66 

50.26 

1.21 


30.002 

15 

73.08 

74.60 

2.08 

53.90 

54.43 

0.98 


30.002 

20 

72.11 

73.07 

1.33 

58.65 

59.03 

0.65 


30.002 

25 

70.92 

71.58 

0.93 

63.79 

64.02 

0.36 


30.002 

30 

69.32 

70.33 

1.46 

68.24 

69.38 

1.67 


30.002 

35 

68.68 

69.48 

1.16 

73.71 

75.10 

1.89 

Average 

30.002 




1.37 



1.13 


34.997 

10 

74.20 

74.83 

0.85 

55.06 

56.43 

2.49 


34.997 

15 

71.98 

73.52 

2.14 

60.35 

61.41 

1.76 


34.997 

20 

70.64 

72.04 

1.98 

65.92 

66.84 

1.40 


34.997 

25 

69.08 

70.61 

2.21 

71.33 

72.68 

1.89 


34.997 

30 

68.44 

69.40 

1.40 

76.57 

78.89 

3,03 


34.997 

35 

67.19 

68.60 

2.10 

82.51 

85.47 

3.59 

Average 

34.997 




1.78 



2.36 


38.274 

10 

73.29 

74.07 

1.06 

59.26 

60.41 

1.94 


38.274 

15 

71.36 

72.79 

2.00 

64.80 

65.90 

1.70 


38.274 

20 

70.17 

71.35 

1.68 

69.82 

71.86 

2.92 


38.274 

25 

69.03 

69.95 

1.33 

76.46 

78.23 

2.31 


38.274 

30 

67.62 

68.78 

1.72 

82.06 

84.99 

3.57 


38.274 

35 

66.54 

68.01 

2.21 

88.14 

92.11 

4.50 

Average 

38.274 




1.67 



2.82 


Table 5. Real and imaginary mean dielectric constants; GW results compared with KS 
model values. 


jj., 



(ii) 


Examining eq(l 1) it is seen that the uncertainty inf' is a result of the uncertainty in A/ and in 

5C . Now the uncertainty in the calibration coefficient comes from two sources: the uncertainty 
in the methanol measurements and the uncertainty in the Gregory and Clarke’s values. This can 
be expressed as 


SC = SC M+ SC rff 


( 12 ) 


where subscript M stands for methanol measurements and ref stands for reference measurements 
of Gregory and Clarke. Here, it is assumed that the methanol measurements made by GW are 
independent of the Gregory and Clarke measurements. Now putting eq( 1 2) into eq(ll) and 
taking the variance of the resulting equations gives 

<t 2 £ . =a 2 ,(SW) + a 2 £ .(C M ) + al(C ref ) (13) 


where 
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(14) 


<7 f ,(SW) = 


_2{C)« 


A/ 


fo 


^■(C„)= 


2(A/)cr c 

/. 


MC*) = 


2{A/)(T ( 

/. 


><•/ 


For a given salinity, the standard deviation of A f is computed separately for each temperature 
value and then the average of these standard deviations is found. This average is denoted by cr A/ . 
The values of Af used in this calculation can be found in Appendix C.2. The quantity <j c is the 
standard deviation of the C coefficients given in yellow in Table 3a and a c is the standard 

deviation of the C constant with respect to an uncertainty of 0.05 in the reference measurements 
of methanol made by Gregory and Clarke. 


Real Dielectric Constant Uncertainty 

Temp 

O't(SW) 

Oe(C) 

CT e (ref) 

a E (total) 

10 

0.25 

0.56 

0.14 

0.62 

15 

0.37 

0.54 

0.13 

0.67 

20 

0.3 

0.53 

0.13 

0.63 

25 

0.44 

0.53 

0.13 

0.7 

30 

0.66 

0.53 

0.12 

0.86 

35 

0.87 

0.51 

0.12 

1.02 

Average 




0.75 


Table 6a. Uncertainty in the real part of the dielectric constant. 


Imaginary Dielectric Constant Uncertainty 

Temp 

<r £ (SW) 

ct £ (C) 

tr E (ref) 

a £ (total) 

10 

0.15 

0.41 

0.1 

0.44 

15 

0.41 

0.45 

0.11 

0.62 

20 

0.17 

0.48 

0.12 

0.53 

25 

0.28 

0.52 

0.13 

0.68 

30 

0.76 

0.58 

0.14 

1.03 

35 

0.77 

0.62 

0.15 

1.06 

Average 




0.73 


Table 6b. Uncertainty in the imaginary part of the dielectric constant. 
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Percent Error in Real Dielectric Constant 

Temp 

<£ f > 

%err(£’) 

%err(C) 

%err(epsref) 

%err(total) 

10 

73.33 

0.34 

0.76 

0.19 

0.85 

15 

71.12 

0.51 

0.76 

0.19 

0.94 

20 

70.2 

0.43 

0.76 

0.19 

0.89 

25 

69.14 

0.64 

0.76 

0.19 

1.01 

30 

69.14 

0.96 

0.77 

0.18 

1.24 

35 

66.41 

1.31 

0.76 

0.19 

1.53 

Average 





1.08 


Table 7a. Total percent error in the real part and its components: variations 
in £ , in the calibration coefficient and in the reference dielectric. 


Percent Error in Imaginary Dielectric Constant 

Temp 

<£”> 

%err(E”) 

%err(C) 

%err(epsref) 

%err(total) 

10 

53.37 

0.27 

0.76 

0.19 

0.83 

15 

59.56 

0.69 

0.76 

0.19 

1.04 

20 

63.42 

0.28 

0.76 

0.19 

0.83 

25 

67.79 

0.41 

0.76 

0.19 

1.01 

30 

76.71 

0.99 

0.76 

0.18 

1.34 

35 

81.83 

0.94 

0.76 

0.19 

1.3 

Average 





1.06 


Table 7b. Total percent error in the imaginary part and its components: 
variations in £* , in the calibration coefficient and in the reference dielectric. 

Following a similar line of reasoning, (^)can be computed from eq(2b) by assuming that C and 
A(l / Q) are independent random variables. The result is: 

(e') = (C)(A{\/Q)) (15) 

and the fluctuations of A(l/g) are given by 

<5A(l/Q) = (c)<5A(l/Q) + {A(l/e))£C ( 16 ) 

Now following the reasoning used for the real part, the variance of the imaginary part of the 
dielectric constant, £'\ can be found. The result is: 
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-o*iSW) + o*(C M ) + o*.(<C«) 


(17) 


where 


ct £ .(5^) = (C><t 4(5 <v(C m )=(a(1/Q))<t C|# <t,.(C^) = (a(1/Q))^ (18) 

For a given salinity, the standard deviation of A(l/Q) is computed separately for each 
temperature value and then the average of these standard deviations is found. This average is 
denoted bycr A0 . The values of A(l/Q) used in this calculation can be found in Appendix C.2. 

An examination of the uncertainties, shown in Tables 6a and 6b, point out that the values tend to 
increase with temperature but the increase is not dramatic. It demonstrates that the uncertainty, 
due to the seawater measurement, increases with temperature. The values in the tables show 
that the contribution of variations in the reference measurements is small. There also does not 
appear to be major differences between the behavior of the uncertainty of the real and imaginary 
parts. The values appearing in Table 6a can be put in terms of percent error by dividing the 
standard deviations by (e and multiplying by 100. The results are shown in Table 7a. The 

value of (s') at each temperature is computed by using eq(10) and averaging A f over all 
measurements with that temperature. The percent error of the imaginary part is computed in a 
similar way by using eq(15) to compute (s")at each temperature. The percent errors for the 

imaginary part are shown in Table 7b. Both tables show that the measurement results, after 
taking all major sources of error into account, vary about 1 % of their mean value. It is observed 
that the calibration measurements of methanol made by GW account for about three quarters of 
the total variation. 

5 Conclusions 

Seawater measurements were made at 30, 35 and 38 psu for a temperature range of 10°C to 35°C 
in 5°C intervals. Several measurements were taken at each salinity and temperature. From these 
measurements, a mean and a standard deviation were calculated. The standard deviation of the 
measurements was lowest at low temperatures (approximately 0.3%) while it was highest at high 
temperatures (approximately 1.0%). In addition to measurement variations, there were variations 
in measurements used to calculate the calibration coefficient. These measurements were taken 
using methanol as a reference liquid at a temperature of 20°C. It was observed that variations in 
the methanol measurements were a substantial contributor to the overall variance of the 
measurement data. Total measurement variation in both the real and imaginary parts due to 
changes in the seawater measurements, in the methanol measurements and in the reference value 
of methanol was less than or equal to approximately 1 .0% for temperatures of 25°C or less. For 
30°C and 35°C, the variation was about 1 .3%. 
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Comparisons of these dielectric measurements with the predictions of the KS model function 
(Klein and Swift 1977) have been discussed. The GW data for the real part of the dielectric 
constant is about 1 .5% below the Klein and Swift data. The difference does not appear to be a 
function of temperature. At present there seems no simple explanation for this difference. The 
imaginary part of the dielectric constant is on the average about 2.0% below the Klein and Swift 
results but the difference tends to increase with both temperature and salinity. This difference 
may be a result of the microwave heating that occurred in the Ho et al data. This heating 
increases with seawater conductivity that, in turn, increases with temperature and salinity. 

In the future measurements attempts will be made to reduce the measurement variations. One 
option is to position the matching networks within the coolant liquid to stabilize their 
temperature. A slightly larger tube radius could be used to increase the frequency difference 
between empty and sample measurements. Assuming that the noise remains the same, this 
would improve the accuracy. The measurements of the calibration liquid can be improved by 
making them at a lower temperature. An examination of the seawater measurements shows that 
the uncertainty diminishes with decreasing temperatures. It is conjectured that the same would 
be true when measuring methanol. In the future, we hope to measure other salinities and to 
measure all salinities at lower temperatures. 
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Appendix A 

Sensitivity of Salinity to Variations in the Model Function 


The following discussion addresses the issue of how accurate the dielectric measurements of 
seawater have to be in order to obtain a given level of uncertainty in the salinity measurements 
from space. Since these dielectric measurements will be used to construct a model function, the 
model function will be as accurate as these measurements . 

To pursue this issue, we start by developing the dependence of the brightness temperature as a 
function of the temperature and the salinity. Consider the relationship between brightness 
temperature and dielectric constant of seawater as given in eq(la) and eq(lb) for the case of a 
radiometer looking directly down, normal to a flat sea surface. By using the equation for the 
model function (eq(lc)) along with eq(la) and eq(lb) we obtain an expression for the brightness 
temperature in terms of the temperature and the salinity as shown below: 

T„ = T H (e(T,S),T) 

Although this relationship has been obtained for a downward pointing radiometer on a flat ocean 
surface, it is valid in much more general conditions. 

It is assumed that the exact model function, £(S,T), is unknown. An approximate model function, 
£i„(S,T), is known. Now if 1'h,> is the observed brightness temperature for seawater at salinity S 0 
and temperature T a , then: 


T Bo = T b (e(T 0 ,S 0 ),T 0 ) (Al) 

Assuming T„ is known exactly and since £,„ is known instead of £ then we naturally understand the 
above equation as: 


T Bo T r (f m (T 0 ,S m ),T 0 ) (A2 ) 

for some value of salinity S,„. 

The question that naturally arises is: How different are S„ and S,„? A quantitative estimate for this 
difference can be obtained from eq(lA) and eq(2A), since both correspond to the same observed 
brightness temperature. Hence one has: 

r,(e(r„s,xr,)=r,(e.<r„s.xr,). <A3> 

Now define the errors Se m and SS m through £ = £ m + d£ m and S m = S 0 + SS m respectively. Note 

that with these definitions 5£,„ corresponds to the modeling error in £ and <5?,,, corresponds to the 
error in retrieved salinity 

Equation (3A) can now be expressed as: 


41 


r,(«.{r.,s.)+&.,r.)=r B (^(r.,s„+«.),r„). 

Expand both sides to a Taylor series in the relevant quantity and assume 3 (Se,„)/dS can be 
neglected. The expanded equation is: 


T B (£ m (T 0 ,S„),T 0 ) + Se' m + ^jjf 


Se I = 


T B (£ m (T 0 ,S 0 ),T 0 ) + 


de' m 


+ 


dT a de. 


B m 


de' 3 5 34 35 


ss„ 


The above expression then yields: 


+ ^-Ss' = 

34 34 


dT B de' m + dT B 3f' 


34 35 34 35 




Hence: 


— V JLSe ' 

op _ ut, m 

m ~ ar B 34 + dr, 34 


34 35 34 35 


Equations (1 ) and (2) can be used to calculate the following derivatives ( £ m »1): 


(A4) 


-^«2Tin 2 Re{4 3/2 } 

3 £ m 

« 27 IT I 2 Im{4 3/2 } 

3 s m 

Inserting the above expressions for the derivatives in eq(4A) yields the error in the retrieved 
salinity as: 




Re |^ 3 ' 2 )&:, + Mf ' 3 ' 2 )&,; 


A worst case error bound is: 




IRe{4 /2 }ll&;,l + IIm{4 3/2 }ll&:[ 

Re{£" 3/2 } ^=- + Im{4 /2 } — 
m 35 '" 35 


(A5) 
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The Klein/Swift model can used to obtain an estimate of the derivatives 3 £' m /dS and d£* n /dS . 


Figure 2 (a) shows the percentage error in retrieved salinity as a function of salinity at various 
temperatures for 0.2 modeling error in both the real and imaginary parts. Figures 2(b and c) show 
that an error in the real part has a larger effect than an error in the imaginary part. Although the 
expression in eq(A5) has been obtained for a special case, the result can to generalized to treat 
oblique observation angles and polarization effects. 


In determining the sample dielectric permittivity using the perturbation technique, uncertainties 
in field quantities established inside the cavity measurement system demand certain 
approximations to be made to these quantities. These approximations, combined with the 
uncertainties in the dimensions of the measurement system, reduce the accuracy of the retrieved 
sample permittivity. Nevertheless the adverse effects of these uncertainties can be mitigated by 
proper calibration. A common calibration procedure involves the use of a reference liquid of 
known dielectric permittivity thus eliminating the need to accurately predict the fields inside the 
cavity as well as the precise knowledge of cavity system dimensions. This calibration procedure 
however has two inherent drawbacks albeit to lesser extent. These drawbacks are (i) the accuracy 
of the calibration is limited by the accuracy of reference liquid dielectric permittivity and (ii) the 
measurement system response is slightly different for the reference liquid and the sample. In this 
section we investigate the accuracy of the perturbation method when used in conjunction with 
the reference liquid calibration procedure. 

Consider a cylindrical cavity having a height l and radius a, with an empty sample tube placed 
along its axis. The cavity is of the transmission type and the sample tube has a very small inner 
radius a, and an outer radius c^. This inner cylindrical region inside the sample tube where the 
liquid sample can be introduced subtends a volume V a . Without the liquid sample, the dielectric 
constant inside the cavity is equal to e o i.e. that of air except for the sample tube whose relative 
dielectric constant is real and is denoted bye,. The relative dielectric profde for the empty cavity 
is denoted by the real profile function £ Kl (r) where r is a radial coordinate. 

The cavity is excited through an aperture A, and the output power is withdrawn by a second 
aperture A 2 . It is assumed that the empty cavity has only wall and end losses and the cavity 
resonates at the radial frequency 0 ) o . In this case it follows from the complex Poynting theorem 
that: 


Appendix B 

Derivation of the Perturbations Equations for the Measurement of 

Dielectric Constants 



(Bl) 


s 
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where {fi„, H a }are the fields inside the empty cavity and S is the inner surface of the cavity with 
unit normal n and ds=nds is the vector differential surface element. Note that the normal 
points out of the cavity. If we denote the total aperture surface by ^4 = .^ uA 2 and the cavity 
walls and ends by S — Ai.e. the cavity inner surface minus the total aperture opening, then eq(l) 
can be expressed as: 

| Steffi,, x Hi }• ds = -J Ste{fi o x H l }• ds (B2) 

S-A A 

This equation expresses the fact that the real power input to the cavity minus the output from the 
cavity through the apertures is dissipated by the lossy walls and ends of the cavity. 

Next consider the cavity with the sample liquid introduced into the sample tube. In this case, 
the relative dielectric profile is the same everywhere except inside the sample tube where the 
relative dielectric constant is now equal to the complex relative dielectric constant of the sample 
i.e.f, = e' +i£ r s . This new relative dielectric profile will be denoted by £ rel (r ) . Hence with this 

notation, £ rel = £ s in V 0 and £ rel — £ rel outside V a . If the fields inside the cavity with the sample 

are now denoted by the pair {E,H}, then again by the complex Poynting theorem it follows 
that: 


J Steffi x H* }• ds + J E ■ E*dv = -J Steffi x H * }• ds 

S-A V a X 


(B3) 


where co is the resonance frequency of the cavity with the sample. Note that f' is the imaginary 
part of the sample relative dielectric constant at co and V n is the volume of the sample. The 

above equation now expresses the fact that the real power input to the cavity through the aperture 
is dissipated by the lossy walls of the cavity and the loss in the sample. 

To a very good approximation we can assume that the wall losses of the cavity do not change 
when the small amount of sample liquid is introduced into the cavity. Formally we can write: 

J 9Ie{£ xH*\ds ~ | ^e{fi 0 x H * }• ds 

S-A S-A 

Then using this approximation in eq(B3) along with eq(B2) leads to the result: 

coeX\ E-E* dv = - J J Steffi x H * }■ ds - J Steffi,, x H * }■ ds (B4) 

Vo U A ) 

implying that the additional real power input to the cavity after the sample is introduced, 
compensates the loss in the sample. 
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Let us now introduce the perturbation fields e — E — E and h = H — H () , the norms of which are 
small compared to the norms of E o and H () . We can then rewrite (B4) as: 

0J£ o £'j E ■ E *dv ~ -J xh* + ex H * }• ds 

V 0 A 

where we have ignored the second order term exh . 

Now let us obtain relations for retrieving the dielectric constant of the sample. For this purpose 
we make use of the vector expression: 

V.(E* o xH+ExH* 0 ) = i(a>-a> 0 )JU 0 H ■ H* o + i ( a>e rt , E-E* a - co o £ rel E-E* o ) (B6) 


where fl n is the magnetic permeability of air. Using the divergence theorem in eq(B6) and 
rearranging yields: 

\[E* e xH+ExH* 0 )-ds + coe/ s \ E ■ E*,dv = ieo£ 0 Ae \e-E* dv 

s V o V t 

-'Atyj (e o e^ E-E*+ fi a H ■ H * ) dv 


where A 0) = Q) n - co and As = e' — 1 . If we now consider the fact that the perturbation fields are 

small, then we can neglect the imaginary parts of all the integrands in the above. Consequently, 
taking the real and imaginary parts of the above yields the following relations respectively: 


and 


J %s(e* x H + E x H* }• ds + co£ 0 £'s\ ^e{£ • E* \dv = 0 


£ AffStef E-E*)dv 

j 1 J 

® j^{£ 0 £ rel E-E* 0+Mll H-H* 0 }dv 

V 


(B7) 


(B8) 


Equation (B8) is a common relation in literature. The real part of the relative dielectric constant 
for the sample can be found by inverting eq(B8) as: 


Ae = 2C 


A co 
O) 


and 


45 


c = 


V 

2 *.Jae{£ •«*}*. 

K, 


The constant C is dependent on the dielectric constant of the sample itself. We can write it as 
the product of a term that depends only the empty cavity features and a term that is also 
dependent on the sample dielectric. For this purpose we make use of the perturbation fields to 
obtain: 


where 


and 


2 %z{EE* 0 }dv "(l + P) 

V 0 

, J k^E, E',+n.H, Hl)dv 
C "” 2 £„•«.**■ 

1 9?e j e ■ E*j ufv 

P=- S_ f * 

I E„ ■ E„ dv 

J o o 


+M,hH’}dv 

v ^ 


(B9) 


(BIO) 


In the above, the constant C Q is a characteristic of the cavity and is independent of the sample. 

We recognize that (W E ) = (1 / A)e n J e rd E 0 • E‘ o dv and (\V M ) = (1 / 4)//„ J // 0 • i/ * rfv are real 

v v 

quantities (assuming negligible loss in the sample tube i.e. e t is real) and respectively correspond 
to time average electrical and magnetic energies stored in the cavity. Since at resonance 

(W F ) = (w u ), we can write J (e„e rel E n • E* + / l a H „ • H* 0 ) dv = S(W E ^. Similarly, the quantity 

V 

(Wto) = C '4Ki E„-E.* dv corresponds to the time average electrical energy stored in the 

K 

empty sample tube. Then the constant C 0 = (W £ )j{W Eo ) is the ratio of the electric energy stored 
inside the sample tube volume to the electric energy stored in the cavity. 

Summarizing the above, we have the relation C = %C 0 . Now we can write Af r as: 
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Ae = 2 xC a ~ 

0) 


(Bll) 


In eq(BI 1)> ;tr = (1 + #)/(l + p)is an ‘energy storage factor’ that reflects the ratio of the electric 
energy stored inside the sample tube volume with the sample present to the total energy stored in 
the cavity relative to the same ratio when the tube is empty. 

Having obtained the real part of the sample dielectric constant we now turn to the imaginary part. 
For this purpose we make use of eq(B7) and put it in a form that involves measured quantities. 
First we rewrite the surface integral in eq(B7) using the perturbation fields as: 

2 1 9te{E 0 x H* \ds + \ Stefe* xh+exH* }■ ds + (0£ o e" s \ ■ K )dv = 0 

s S v 0 

Realizing from eq(Bl) that the first term on the left disappears and using eq(B5) in the above, we 
recognize that 

1 • E* \dv * J ^{e E*)dv 

V V 

v o v o 

to the first order. Then using eq(B4) we can write 

<D£ 0 £ s\ \dv = -j J 5Ke{jE' x H * }■ ds - J x H* }• ds i (B 1 2) 

K U A \ 

Note that we can express the integrals on the right hand side of the above in terms of the Q- 
factor. The Q-factor is defined as: 


_ ^ Time Average Energy Stored 

Time Average Energy Loss / second 


Using this definition, we can write the second term on the right hand side of eq(B12) with the 
help of eq(B2) as: 


~ J x }■ ds =£rlb £ ”' E - • E • + v-"" ■ "* )* 

z a v 



(B13) 


Similarly with the help of eq(B3), the first term on the right hand side of eq(BI2) can be written 
as: 


Z A 






(B 14) 
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Note that 


l Mefo, }E-E*+fii 0 H-H*)dv = l (e o e rel E -E* +/x 0 H ■H*)dv + e g Aej v EE*dv 


Using the perturbation fields and eq(B9) and eq(B10), to the first order, the integrals on the right 
hand side of this expression can be written as: 


J E ■ E*dv « (1 + 2p)j E 0 ■ E*dv 


and 


\k 


e rel E-E' 


+ m„H • H * )dv - (1 + 2 q)\ (e o £ rel E 0 ■ E* 0 +fi 0 H 0 ‘ H* a )dv 

V V 

Using these results along with eq(Bl 1), we can write eq(B14) as: 


-U*t\pxH*\ds = 2^{W E )[\ + 2q + (\ + 2p)x 
^ a y \ 


A CO 

(0 J 


Since A col 0) < 10 4 , we can ignore the third term in the brackets on the right hand side: 


1 X H* }• « 2— {W h: )(l + 2 q) 

A Q 


(B15) 


Using eq(B13) and eq(B15) in eq(B12) we obtain: 


'f9?efc-<}rfv»4 

v j___0 

+ 

2 q 

V 0 

yQ Q 0 j 


Q_ 




In the above we have once again neglected a term of order A 0)1(0 . Recognizing that 
£„ £ E ■ E*dv = 4(1 + p)(W Eo ) , we obtain: 


£. -■ 


c 0 

v i 

_ n 

+ 2 £ 

0 + p) 

[q 

Qo j 

q\ 


If we add and subtract q/Q 0 to the right hand side, we finally obtain 

, i V 


< “ %Co 


f \ 1 ^ 


Q Qo) 0 + «) 


Q + Q„ 


(B 16 ) 
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This equation is more accurate than the commonly recognized form which can be obtained by 
ignoring the second term in eq(B16) as: 


where A(l/0 = {MQ-MQ ,,) . 


: xc,A 




Q 


\x-j 


(B 17) 


ERROR ANALYSIS 

The constant C g is very sensitive to the inner volume of the sample tube and is thus hard to 
predict accurately. However by using a reference liquid with known dielectric constant 
£ r = e' + ie " , the constant C may be eliminated as follows: 


~ 1 _ X Ar» 

e’r - I Xr ^°>r 


(B18) 


where % r = (1 + q r ) /(1 + p r ) and p r and q r are obtained from eq(B3) and eq(B4) respectively 
when the reference liquid is used, and Ary,. is the frequency shift corresponding to the reference 
liquid. Similarly assuming we can suffice with the approximate expression given by eq(B17), we 
can write for e' : 


< = x A(l/g) 

< ^ A(1 1 Q r ) 


(B19) 


Where A(1 / <2 r ) = (1/Q,. — \l Q„) and Q r is the Q-factor of the cavity when the reference liquid is 
introduced. In practice, instead of estimating %and X r * the approximations 


e' s -1 _ o) r A a 
£ r -1 0) A Q) r 

< ^ A(1 / Q) 

< A(1 !Q r ) 


(B20) 


(B21) 


are used. This can be justified on the basis of the fact that \p\ « \,\q\ « l,|p r | « 1 and \q r \ « 1 . 
The above approximation yields better results if the dielectric constant of the reference liquid is 
close to that of the sample liquid. The dimensionless quantity \x! X r — 1| is a direct measure of 

the fractional error in sample dielectric permittivity incurred in the perturbation approximations 
of eq(B20) and eq(B21). In the following we estimate this error and present its dependence on 
the sample tube diameter. 
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In order to estimate the error we have to determine the constants p,q,p r and q r . For this 
purpose we begin with the expression of the fields inside the cavity given as: 


and 


E = 


H = 


A Ao( ky[e^p)z 

^ Aj 2*^0 ^ o 2^0 [kijsi 

[A o3 J 0 {kp) + B o3 N 0 {kp)]z 

n;' k 2 i, (k^p)+ b o 2 n , (kjppp 

TJ- l [A oi J ] (kp)+B o3 N l (kpM 


0 < p < a, 

a, < p<a 2 
a 2 < p <a 3 


0< p<a t 
a x < p <a 2 
a 2 <p<a 3 


where q rel = \f . For the empty cavity one has k = k 0 = co 0 A„P 0 and £rel = 1 . Similar 
expressions are valid when the sample liquid is inserted into the sample tube. In this case, 
k = k n + Act)/(c-yje^) , £ rel =e x and the coefficients A ni and B m are replaced by A si and B si 

respectively for i = {1,2,3}. In the case of the reference liquid, k =k a +A o> r /{cyf£^),£ rc i = £ r , 
and the coefficients A nj and B oi are replaced by A n and B ri respectively for (' = {1,2,3}. The 
coefficients are found by applying the boundary conditions i.e. the continuity of E and H fields 
at p = a x ,p = a 2 and p - a 3 (note that both field have only tangential components to the 

boundaries). It is assumed that the walls and end plates are perfectly conducting. This leads to 
the following relations 


A , 


A . 


A 


A tl 


A„. 


Fa ,1 


= [rj 

01 

9 

St 

= [r,l 


and 

ri 

= [r„l 

rl 

A. 


A. 


A. 


A. 


A. 


A. 


where [rj , [rj and [r ,] are the transition matrices for the empty cavity, cavity with sample 

liquid and cavity with reference liquid respectively, relating the coefficients of the V th layer to 
those of the first layer. The generic transition matrix is given by 


where 


[r,]=nW 

P=1 
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M 


7lX L ’ P 

JLA - p+\,p 


n p N 0 (X p+Up )J l (X P ' P )-N l (X p+lp )J 0 (X P ' P ) 

L-«/o(^, + , P )-/,(^ P .p) + A(^, +I . P )-/o(^,, P ) - n pJo(X p+i p ) N { (X p p ) + J, (^p+i.p ) N 0 (X p p )- 


with X p? = k 3 Je^a q , n p =^e^Je~ x and£ p is the relative dielectric permittivity of the p’th 
layer. 

Using the field expressions, we can now obtain explicit expressions for the constants p,q , p r and 
q r . However, in order to obtain simple expressions we make several approximations. Since a , is 
so small that k 3 fe^a l « 1 , we can write J 0 {{k g + A(oj(cy[e^)}yf£~ s p] = J 0 (k (l p) ~ 1 . Hence, in 

the region 0</7<a,we have e = A sl J 0 (k s p)z- A ol J 0 (k o p)z ~ (A s] - A oX )z . Thus from eq(B9) 

_ g„ (9te{ A v] }/ 4, - 1)V, _ 9te{ A„ } ] 

e. Vi A,> 


and consequently 1 + p ~ 9te{A s] }/A o] for the sample liquid. Similarly, for the reference liquid 
we obtain 1 + p r ~ 9te{A rl }/A„, . For the quantity q given in eq(B10), we first note that 
a { ,a 2 « a 3 , hence eq(B10) can be approximated by 


J9 le{ £n e rel e-E*+ Mo h-H*}dv 



\{e 0 £ rel E 0 .E* 0 +M 0 H 0 -H* 0 )d V 


Since A COa 3 /c « 1, we can write J n ([/:„ +A col c\p)~ J „{k n p) in the region a 2 < p <a 3 . This 


approximation leads to 
pertinent region. Thus 


e “ (^*3 ~ A o3 )J 0 (k o p)z and h~T]; l (A s3 - A oZ )N x (k oP )<j> 
(9te{ A j3 } - A n3 )/ 


q 


■v3 J / x n?> J 

A** 


= 9te{A j3 }/A o3 - 1 


in the 


where I = f ’ p(e 0 J l (k 0 p) + r/ 0 i p 0 J?(k 0 p))dp . With this final expression we obtain the result 

Ja 2 

• • 

\ + q ~ 9?e{ A t3 }/A o3 . Note that A o3 /A nl = [r„ 3 ],,=/„ (real quantity) and A s 3 /A sl = [T s3 ]„ = y s . 
Hence with the help these relations, we finally obtain 2'~(^ e {A s3 }/9te{A sl })(A fll /A fj3 ) 
~ 9tc{ A 0 /A sl }(A o1 / A o3 ) = 9te{ y s }/ y o . In this final expression, the second approximation is based 
on the premise that in V 3 the fields are negligibly perturbed when the sample is introduced so 
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that these fields are almost real. Similarly, defining y r =A rl l A rl ~ [r r3 ] M for the reference liquid 
we obtain x r ~ ^e{ y, > / y o . With this final, relation the fractional error reduces to 
\x'X r -l|»|9te{ft}/9le{r r }-l|. 
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Appendix C 
Measurement Data 


This appendix contains three subsections. In the first subsection, the network configuration 
(shown in Figure 7) is given for the methanol and seawater measurements in four separate tables. 
In subsection C.2 the seawater measurement data is given is three tables. Finally, in subsection 
C.3, a table of the frequency drift over the time span of an experiment is given for several 
representative cases. 

C.l Configuration Data 


Methanol Measurement Configuration 

Date 

Exp# 

Methanol 

Tube 

Temp 

Isolator 

Attenuator 

Amplifier 

01/26/07 

1 

New 

#12 

20 °C 

Yes 

No 

Yes 

01/26/07 

2 

New 

#12 

20 °C 

Yes 

No 

Yes 

01/26/07 

3 

Old 

#12 

20 °C 

Yes 

No 

Yes 

01/27/07 

1 

New 

#12 

20 °C 

Yes 

No 

Yes 

02/01/08 

1 

New 

#12 

20 °C 

Yes 

No 

Yes 

02/04/08 

2 

New 

#12 

20 °C 

Yes 

No 

Yes 

04/1 0/08 

1 

New 

#12 

20 °C 

Yes 

Yes 

No 

04/19/08 

i 

New 

#12 

20 °C 

Yes 

Yes 

No 

04/24/08 

i 

New 

#12 

20 °C 

Yes 

Yes 

No 

04/24/08 

2 

New 

#12 

20 °C 

Yes 

Yes 

No 

04/25/08 

1 

New 

#12 

20 °C 

Yes 

Yes 

No 

04/29/08 

1 

New 

#12 

20 °C 

Yes 

Yes 

No 

04/30/08 

i 

New 

#12 

20 °C 

Yes 

Yes 

No 

04/30/08 

2 

New 

#12 

20°C 

Yes 

Yes 

No 

05/13/08 

1 

New 

#12 

20 °C 

Yes 

Yes 

No 


Table C. 1-1 . Configuration for Methanol Measurements 
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Seawater Measurement Configuration (30.002 psu) 

Date 

Tube 

Temp 

Exp# 

Isolators 

Attenuators 

Amplifier 

Freq Stabilizer 

12/14/07 

#12 

10 

1 

Yes 

No 

Yes 

Yes 

12/28/07 

#12 

10 

1 

Yes 

No 

Yes 

Yes 

12/28/07 

#12 

10 

2 

Yes 

No 

Yes 

Yes 

02/08/08 

#12 

10 

i 

Yes 

No 

Yes 

Yes 

12/29/07 

#12 

15 

1 

Yes 

No 

Yes 

Yes 

12/29/07 

#12 

15 

2 

Yes 

No 

Yes 

Yes 

02/07/08 

#12 

15 

1 

Yes 

No 

Yes 

Yes 

12/09/07 

#12 

20 

i 

Yes 

No 

Yes 

Yes 

12/12/07 

#12 

20 

i 

Yes 

No 

Yes 

Yes 

12/29/07 

#12 

20 

i 

Yes 

No 

Yes 

Yes 

12/30/07 

#12 

20 

i 

Yes 

No 

Yes 

Yes 

12/30/07 

#12 

20 

2 

Yes 

No 

Yes 

Yes 

02/07/08 

#12 

20 

i 

Yes 

No 

Yes 

Yes 

1 2/26/07 

#12 

25 

1 

Yes 

No 

Yes 

Yes 

12/28/07 

#12 

25 

i 

Yes 

No 

Yes 

Yes 

12/30/07 

#12 

25 

1 

Yes 

No 

Yes 

Yes 

02/05/08 

#12 

25 

1 

Yes 

No 

Yes 

Yes 

02/05/08 

#12 

25 

2 

Yes 

No 

Yes 

Yes 

02/06/08 

#12 

25 

1 

Yes 

No 

Yes 

Yes 

12/26/07 

#12 

30 

1 

Yes 

No 

Yes 

Yes 

12/27/07 

#12 

30 

1 

Yes 

No 

Yes 

Yes 

02/10/08 

#12 

30 

1 

Yes 

No 

Yes 

Yes 

12/15/07 

#12 

35 

1 

Yes 

No 

Yes 

Yes 

12/27/07 

#12 

35 

1 

Yes 

No 

Yes 

Yes 

12/27/07 

#12 

35 

2 

Yes 

No 

Yes 

Yes 

02/1 0/08 

#12 

35 

i 

Yes 

No 

Yes 

Yes 


Table C.I-2. Network Configuration for Seawater Measurements with Salinity 30.002 psu 


Seawater Measurement Configurate 

on (34.997 psu) 

Date 

Tube 

Temp 

Exp # 

Isolators 

Attenuators 

Amplifier 

Freq Stabilizer 

05/20/08 

#12 

10 

1 

Yes 

Yes 

No 

Yes 

05/30/08 

#12 

10 

1 

Yes 

Yes 

No 

Yes 

05/20/08 

#12 

15 

1 

Yes 

Yes 

No 

Yes 

05/30/08 

#12 

15 

1 

Yes 

Yes 

No 

Yes 

05/1 9/08 

#12 

20 

1 

Yes 

Yes 

No 

Yes 

05/29/08 

#12 

20 

1 

Yes 

Yes 

No 

Yes 

05/19/08 

#12 

25 

1 

Yes 

Yes 

No 

Yes 

05/28/09 

#12 

25 

1 

Yes 

Yes 

No 

Yes 

05/27/09 

#12 

30 

1 

Yes 

Yes 

No 

Yes 

05/28/08 

#12 

30 

1 

Yes 

Yes 

No 

Yes 

06/02/08 

#12 

30 

1 

Yes 

Yes 

No 

Yes 

06/04/08 

#12 

30 

1 

Yes 

Yes 

No 

Yes 

05/27/08 

#12 

35 

i 

Yes 

Yes 

No 

Yes 

05/29/08 

#12 

35 

1 

Yes 

Yes 

No 

Yes 


Table C.l-3. Network Configuration for Seawater Measurements with Salinity 34.997 psu 
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Seawater Measurement Configuration (38.274 psu) 

Date 

Tube 

Temp 

Exp # 

Isolators 

Attenuators 

Amplifier 

Freq Stabilizer 

03/05/08 

#12 

10 

1 

Yes 

No 

No 

Yes 

05/02/08 

#12 

10 

1 

Yes 

Yes 

No 

Yes 

03/04/08 

#12 

15 

1 

Yes 

No 

No 

Yes 

03/05/08 

#12 

15 

1 

Yes 

No 

No 

Yes 

05/02/08 

#12 

15 

1 

Yes 

Yes 

No 

Yes 

03/05/08 

#12 

20 

1 

Yes 

No 

No 

Yes 

05/01/08 

#12 

20 

1 

Yes 

Yes 

No 

Yes 

05/08/08 

#12 

20 

1 

Yes 

Yes 

No 

Yes 

03/03/08 

#12 

25 

1 

Yes 

No 

No 

Yes 

05/01/08 

#12 

25 

1 

Yes 

Yes 

No 

Yes 

03/03/08 

#12 

30 

1 

Yes 

No 

No 

Yes 

05/05/08 

#12 

30 

1 

Yes 

Yes 

No 

Yes 

05/06/08 

#12 

30 

1 

Yes 

Yes 

No 

Yes 

05/14/09 

#12 

30 

1 

Yes 

Yes 

No 

Yes 

05/23/08 

#12 

30 

1 

Yes 

Yes 

No 

Yes 

03/04/08 

#12 

35 

1 

Yes 

No 

No 

Yes 

05/05/08 

#12 

35 

1 

Yes 

Yes 

No 

Yes 

05/06/08 

#12 

35 

1 

Yes 

Yes 

No 

Yes 

05/15/08 

#12 

35 

1 

Yes 

Yes 

No 

Yes 

05/23/08 

#12 

35 

1 

Yes 

Yes 

No 

Yes 


Table C.l-4. Network Configuration for Seawater Measurements with Salinity 38.274 psu 
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C.2 Seawater Dielectric Constant Data 
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Table C.2-3. Seawater Dielectric Measurement for Salinity at 38274 psu 


C.3 Frequency Stability 


In this subsection of Appendix C, the long term frequency stability of the system is discussed. To 
help with the discussion, a frequency stability table has been constructed. In this table the 
difference between the start resonant frequency and the stop resonant frequency is recorded for 
1 8 different experiments. The time difference is also recorded. Here it is seen that the average 
frequency drift over one and a half hours is approximately 1 KHz. Since the cavity is held at a 
stable temperature (±0.02°C), it is believed the reasons for this drift are changes is the room 
temperature. There is, however, no observable correlation between changes in the room 
temperature and changes in the cavity resonant frequency. This is because the temperature 
changes in the room have an integrated effect on the network external to the liquid bath. 
Although these resonant frequency changes occur, they do not have a noticeable effect on the 
measurements since they occur slowly. They do however, make it necessary to do the 
measurements one at a time rather than automating the system. This has seriously slowed down 
the amount of data we could collect. One proposed solution has been to put the isolators and 
attenuators in the liquid bath. Their temperatures would be held constant and thus the effects of 
external temperature changes would be much reduced. 


Frequency Stabilil 

ty (Seawater) 

Date 

psu 

Temp 

(°C) 

Att. 

Start Freq.(Hz) 

Start Time 

End Freq.(Hz) 

End Time 

Af(Hz) 

At 

12/28/07 

30.002 

10 

N 

1,413,492,578 

10:49 

1,413,491,691 

13:03 

887 

2:13 

12/29/07 

30.002 

15 

N 

1,413,492,578 

4:24 

1,413,357,723 

7:22 

2729 

2:58 

02/07/08 

30.002 

20 

N 

1,413,492,578 

14:45 

1,413,249,334 

18:05 

1523 

3:19 

12/30/07 

30.002 

25 

N 

1,413,492,578 

18:48 

1,413,092,123 

20:35 

326 

1:46 

12/26/07 

30.002 

30 

N 

1,413,492,578 

9:38 

1,412,961,884 

12:25 

38 

2:47 

12/27/07 

30.002 

35 

N 

1,413,492,578 

4:25 

1,412,832,612 

6:55 

2283 

2:29 

05/30/08 

34.997 

10 

Y 

1,413,492,578 

16:44 

1,413,450,930 

18:55 

1209 

2:11 

05/20/08 

34.997 

15 

Y 

1,413,492,578 

10:58 

1,413,339,361 

12:33 

485 

1:35 

05/29/08 

34.997 

20 

Y 

1,413,492,578 

12:57 

1,412,781,762 

14:36 

2521 

1:38 

05/19/08 

34.997 

25 

Y 

1,413,492,578 

11:49 

1,413,076,114 

12:57 

157 

1:07 

05/27/08 

34.997 

30 

Y 

1,413,492,578 

12:57 

1,412,920,989 

14:10 

55 

1:13 

05/27/08 

34.997 

35 

Y 

1,413,492,578 

15:09 

1,412,789,412 

16:36 

2234 

1:27 

05/02/08 

38.274 

10 

Y 

1,413,492,578 

16:57 

1,413,473,408 

18:44 

534 

1:46 

03/04/08 

38.274 

15 

N 

1,413,492,578 

10:56 

1,413,357,029 

12:05 

1491 

1:09 

05/01/08 

38.274 

20 

Y 

1,413,492,578 

10:52 

1,413,221,754 

12:36 

355 

1:43 

05/01/08 

38.274 

25 

Y 

1,413,492,578 

15:05 

1,413,098,226 

16:35 

733 

1:30 

05/06/08 

38.274 

30 

Y 

1,413,492,578 

11:20 

1,412,949,791 

12:56 

40 

1:35 

05/06/08 

38.274 

35 

Y 

1,413,492,578 

14:49 

1,412,818,111 

16:31 

2528 

1:42 

Average 








1118 

1:54 

Standard 

Deviation 








970 

0:38 


Table C.3-1 Frequency Stability Table. Start and end times occur on the same day. 
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Appendix D 

Loaded and Unloaded Q 


In this appendix, the subject of loaded and unloaded Q will be discussed. The formulas for the 
imaginary dielectric constant given in eq(2a) assume that the Q used is the unloaded Q. It will be 
denoted by Qul- The Q measurement taken by the NA is the loaded Q that will be denoted by 
Q l . For transmission cavities the two are related by 

Qul = Ql( 1+ P\ + A) (Di) 

where and f3 2 are the input and output coupling coefficients. For a discussion of these 
quantities, please refer to the “Measurement of Q” by M. Sucher (1963) or Kajfez et al. (1999). 
If the cavity is viewed as a reciprocal two port, its scattering parameters are S n ,S 21 and S n . The 

scattering parameter S 2I is directly related to the coupling coefficients. At the resonant frequency 
of the cavity, the relationship is given by 


S 


21 


1 + P\ + A 


(D2) 


If it is assumed that the cavity coupling coefficients are the same (approximately true for our 
cavity) then /? = ft = /? 2 and eq(D2) can be solved for /? , thus 

2jB = 7ZV~ (D3) 

1 i 2I 

Equations (Dl) and (D3) can be applied to our Q measurements to determine Qul from Qi.. First, 
it is noted that the magnitude of S 21 is measured by the NA between terminals a and b of Figure 
7. This is not the S 21 of the cavity, but knowing the network configuration given in Appendix 
C. 1, the cavity S 21 can be calculated. When only isolators are used, 1 dB is subtracted from the 
S 21 computed by the NA while when both isolators and attenuators are used , 1 1 dB is subtracted. 

Some results of these calculations are presented in Table D-l for the empty cavity Q and for the 
Q when the cavity has a sample. The columns marked (S ) dB and (S 2I ) dB are the cavity values 

of S 21 in dB for the empty cavity and for the cavity with a sample in it respectively. An 
examination of the table shows that the loaded and unloaded Q’s can differ by as much 300 to 
500. However, when A(l/Q)is computed for the loaded and unloaded cases, there is very little 
difference as demonstrated by the %Error column. The errors are under 0.1 % in general and thus 
there is no need to compute the unloaded Q when calculating the imaginary part of the sample 
dielectric constant from eq(2a). 
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LOADED and UNLOADED Q 

Date 

psu/Meth. 

Temp (C) 

Att 

(Lnw)cJB 

Ql E 

(S21 E )dB 

Qul E 

Ql S 

(S21 S )dB 

Qul S 

A(1/Q L ) 

A (1/Qui_) 

%Error 

02/08/08 

30.002 

10 

N 

1 

7595 

-26.15 

8039 

5130 

-29.54 

5329 

6.327E-05 

6.325E-05 

0.02% 

12/29/07 

30.002 

15 

N 

1 

7504 

-24.99 

8010 

4980 

-28.66 

5195 

6.754E-05 

6.765E-05 

0.16% 

12/30/07 

30.002 

20 

N 

1 

7577 

-27.35 

7960 

4851 

-31.24 

5005 

7.416E-05 

7.41 7E-05 

0.01% 

12/30/07 

30.002 

25 

N 

1 

7733 

-29.65 

8030 

4731 

-33.84 

4841 

8.206E-05 

8.201 E-05 

0.05% 

12/27/07 

30.002 

30 

N 

1 

7693 

-28.51 

8031 

4603 

-33.03 

4721 

8.726E-05 

8.729E-05 

0.04% 

12/28/07 

30.002 

35 

N 

1 

7647 

-26.02 

8101 

4462 

-30.73 

4613 

9.334E-05 

9.337E-05 

0.02% 

05/20/08 

34.997 

10 

Y 

11 

7840 

-30.14 

8813 

5070 

-33.94 

5459 

6.969E-05 

6.972E-05 

0.05% 

05/20/08 

34.997 

15 

Y 

11 

7837 

-30.33 

8786 

4906 

-34.41 

5261 

7.623E-05 

7.625E-05 

0.02% 

05/29/08 

34.997 

20 

Y 

11 

7830 

-31.29 

8669 

4736 

-35.63 

5031 

8.343E-05 

8.340E-05 

0.04% 

05/19/08 

34.997 

25 

Y 

11 

7834 

-31.16 

8687 

4591 

-35.80 

4871 

9.01 7E-05 

9.017E-05 

0.00% 

05/28/08 

34.997 

30 

Y 

11 

7821 

-31.54 

8632 

4434 

-36.40 

4686 

9.767E-05 

9.757E-05 

0.10% 

05/27/08 

34.997 

35 

Y 

11 

7820 

-32.62 

8528 

4308 

-37.90 

4512 

1.042E-04 

1.044E-04 

0.12% 

05/02/08 

38.274 

10 

Y 

11 

7382 

-29.98 

8317 

4743 

-33.84 

5112 

7.537E-05 

7. 539 E-05 

0.02% 

05/02/08 

38.274 

15 

Y 

11 

7433 

-29.80 

8397 

4628 

-33.93 

4984 

8.154E-05 

8.156E-05 

0.02% 

05/01/08 

38.274 

20 

Y 

11 

7534 

-29.80 

8512 

4522 

-34.28 

4855 

8.841 E-05 

8.850E-05 

0.10% 

03/05/08 

38.274 

10 

N 

1 

7684 

-34.28 

7854 

4875 

-38.21 

4943 

7.499E-05 

7.498E-05 

0.01% 

03/05/08 

38.274 

20 

N 

1 

7702 

-34.25 

7873 

4576 

-38.79 

4636 

8.870E-05 

8.870E-05 

0.01% 

03/03/08 

38.274 

25 

N 

1 

7596 

-31.10 

7841 

4378 

-35.82 

4459 

9.677E-05 

9.674E-05 

0.03% 

03/03/08 

38.274 

30 

N 

1 

7542 

-30.52 

7803 

4195 

-35.60 

4275 

1.058E-04 

1.058E-04 

0.01% 

01/26/08 

Meth. 

20 

N 

1 

7741 

-25.25 

8246 

6979 

-26.17 

7386 

1.409E-05 

1.412E-05 

0.16% 

01/26/08 

Meth. 

20 

N 

1 

7752 

-25.41 

8249 

6981 

-26.34 

7381 

1.424E-05 

1.426E-05 

0.11% 

01/27/08 

Meth. 

20 

N 

1 

7779 

-26.14 

8234 

7004 

-27.07 

7370 

1.423E-05 

1.424E-05 

0.10% 

04/30/08 

Meth. 

20 

Y 

11 

7569 

-29.74 

7857 

6838 

-30.60 

7073 

1.412E-05 

1.41 IE-05 

0.07% 

04/30/08 

Meth. 

20 

Y 

11 

7573 

-29.76 

7860 

6833 

-30.64 

7066 

1.430E-05 

1 .429E-05 

0.04% 

05/13/08 

Meth. 

20 

Y 

11 

7846 

-30.50 

8118 

7068 

-31.40 

7288 

1.404E-05 

1.404E-05 

0.03% 


Table D-l Loaded and Unloaded Q values 
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